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1.0  INTRODUCTION 


1.1  Significance  of  the  Problem 

Advanced  materials  for  aerospace,  structural,  power  and  propulsion  applications  offer 
significant  advantages  in  terms  of  efficiency  and  cost.  This  has  been  the  reason  for  the  ever 
increasing  use  of  advanced  composite  materials  in  recent  years.  For  example,  most  new  or 
recently  manufactured  aircraft  include  structural  pans  made  of  composite  materials.  Although 
practical  usage  of  composite  materials  in  the  modem  age  originated  in  aircraft  and  aerospace 
industries,  their  advantages  are  now  recognized  by  all  major  industries.  A  widespread  and 
efficient  application  of  composite  materials  requires  detailed  and  reliable  knowledge  ci  th?’r 
physical  propenies  and,  in  turn,  of  their  behavior  under  applied  loads.  There  are  a  number  of 
important  technical  problems  associated  with  the  use  of  composite  materials.  One  such  prob¬ 
lem  is  the  effect  of  discontinuities  (holes  and  notches)  on  the  strength  of  composite  laminates. 
This  issue  is  critical  for  the  determination  of  the  load  bearing  capacity  of  composite  laminates; 
which  is  directly  applicable  to  the  design  of  composite  panels  and  the  location  of  fastener 
holes.  Indeed,  the  manufacture  and  repair  of  advanced  composite  stmctures  have  serious 
problems  connected  with  the  placement  of  fastener  holes.  This  is  especially  relevant  to  com¬ 
posite  panel  repair,  both  in  the  field  and  at  the  repair  facility.  At  the  present  time  all  depots 
are  confronted  with  these  problems.  The  lack  of  appropriate  data  has  resulted  in  new  and  in- 
service  designs  which  are  often  unnecessarily  conservative  and  expensive  (both  in  cost  and 
tum-around-time).  Another  related  problem  concerning  composite  materials  is  the  issue  of 
interlaminar  response  of  composite  materials  which  is  directly  related  to  delamination  and 
edge  effects  in  composites.  In  recent  time,  delamination  has  become  the  most  feared  failure 
mode  in  laminated  composite  structures.  It  can  exhibit  unstable  crack  growth,  and  while 
delamination  failure  itself  is  not  usually  a  catastrophic  event,  it  can  perpetrate  such  a  condi¬ 
tion  due  to  its  weakening  influence  on  a  component  in  its  resistance  to  subsequent  failure 
modes.  What  had  begun  in  1970  as  somewhat  of  an  academic  curiosity  turned  into  a  beehive 
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of  research  activity  in  recent  years.  Study  of  delamination  is  now  one  of  the  most  prominent 
topics  in  composite  mechanics  research.  Another  related  issue  in  the  engineering  application 
of  composite  materials  is  the  use  of  numerical  methods  and,  in  particular,  the  finite  element 
method  in  linear/nonlinear  modeling  of  composite  laminates.  The  various  composite  shell  and 
composite  solid  elements  that  are  available  today  are  not  adequate  for  advanced  applications. 
These  elements  are  formulated  using  one  or  another  form  of  classical  shell  theories  and 
although  they  provide  acceptable  results  in  simple  loading  conditions,  but  they  are  not  capable 
of  predicting  accurate  structural  response  for  extreme  in-use  loading  conditions  of  aerospace 
systems.  One  reason  for  this  deficiency  has  been  the  lack  of  powerful  theory  for  composite 
laminates.  Yet  another  issue  in  composite  laminates  is  the  development  of  constitutive  rela¬ 
tions  which  adequately  represent  the  effect  of  the  constituents.  In  general,  our  knowledge  of 
the  thermo-mechanical  behavior  of  composite  materials  and  their  constitutive  relations  has 
lagged  behind  advances  in  such  other  related  areas  as  increasingly  sophisticated  computers 
and  computational  methods  for  solution  of  complex  problems.  Proper  understanding  and 
characterization  of  available  and  new  composite  laminates  and  implementation  of  their  consti¬ 
tutive  models  in  modem  (computational)  solution  procedures  is  a  vital  component  for  safe, 
economical  and  competitive  design  in  aerospace  applications.  It  is  believed  that  for  a  reliable 
structural  analysis,  the  constitutive  models  of  composite  materials  should  be  based  on  sound 
thermo-mechanical  principles.  All  the  foregoing  problems  share  one  common  deficiency, 
namely,  the  lack  of  an  adequate  and  sound  theory  applicable  to  composite  laminates,  a  theory 
predicated  upon  solid  physical  principles  and  sound  mathematical  analysis  that  can  account  for 
the  effects  of  micro- structure,  material  nonlinearity,  geometric  nonlinearity,  interlaminar 
stresses,  and  complex  geometry. 

In  phase  I  of  this  research  project  the  feasibility  of  the  development  of  such  a  theory  was 
studied.  The  present  work  is  a  continuation  of  the  phase  I  effort  to  develop  a  complete 
thermo-mechanical  theory  for  composite  laminates  that  exhibits  the  following  characteristics; 
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a)  It  accounts  for  the  effects  of  micro- structure, 

b)  It  accounts  for  the  effects  of  geometric  nonlinearity. 

c)  It  accounts  for  the  effects  of  curvature. 

d)  It  accounts  for  the  effects  of  interlaminar  stresses.  The  three  components  of  the 

interlaminar  stresses  are  included  in  the  theory  and  can  be  determined  numerically. 

e)  It  has  a  continuum  character. 

f)  It  is  applicable  to  both  static  and  dynamic  problems. 

Due  to  the  use  of  Cosserat  surface  theory,  we  have  called  the  proposed  theory  Cosserat  com¬ 
posite  theory. 

In  the  past,  several  theories  have  been  proposed  for  the  modeling  of  multilayered  plates 
and  shells.  Noor  and  Burton  (1990)  assessed  the  computational  models  for  multilayered  com¬ 
posite  shells.  They  presented  a  list  of  400  related  publications  and  they  identified  the  follow¬ 
ing  four  general  approaches  for  constructing  two-dimensional  theories  for  multilayered  shells: 

1.  method  of  hypotheses. 

2.  method  of  expansion. 

3.  asymptotic  integration  technique. 

4.  iterative  methods  and  methods  of  successive  corrections. 

• 

Some  of  these  approaches  were  reveiwed  in  the  phase  I  of  this  research  project.  Although 
each  of  these  theories  has  its  own  advantages,  we  v/ould  like  to  emphasize  that,  as  a  general 
assessment,  none  of  these  theories  possesses  the  features  of  Cosserat  composite  theory  collec¬ 
tively.  In  particular,  the  inteilaminar  stresses  were  not  addressed  in  most  of  these  theories  and 
their  applications  were  limited  to  situations  where  the  composite  laminate  is  made  of  a  few 
number  of  layers.  These  limitations  are  not  present  in  the  Cosserat  composite  theory. 


-4- 


The  Cosserat  composite  theory,  because  of  its  continuum  character,  is  best  suited  for 
modeling  thick  composite  laminates  composed  of  several  (order  of  hundreds)  thin  plies. 
Furthermore,  the  interlaminar  stresses,  which  are  responsible  for  delamination  failures,  are 
incorporated  into  the  formulation  of  the  theory  in  a  natural  and  consistent  manner  and  without 
any  ad  hoc  assumptions.  Considering  these  technical  issues,  the  objectives  of  our  research 
project  were  developed.  In  the  following  a  detailed  description  of  these  objectives  is 
presented. 

1.2  Objectives  of  the  Present  Research  Project 

Some  technical  issues  associated  with  the  application  of  composite  laminates  were 
addressed  above.  The  present  research  effort  is  directed  toward  studying  these  issues.  In  this 
regard  the  following  objectives  have  been  identified  to  be  addressed  in  the  project: 

1.  Development  of  strain  measures  based  on  the  kinematics  of  micro-  and  macro- 

( 

structures  of  composite  laminates  and  relating  these  strain  measures  to  the  compo¬ 
site  stress  and  composite  stress  couple  tensors.  In  particular,  development  of  consti¬ 
tutive  relations  for  elastic  behavior  of  composite  laminates  in  the  context  of 
infinitesimal  deformations  that  can  be  applicable  to  any  loading  condition  including 
pure  bending  and  pure  extension  is  a  part  of  this  task. 

2.  Further  development  of  the  Cosserat  composite  theory  to  a  thermo-mechanical 
theory  for  multi-constituent,  anisotropic  laminates.  The  theory  outlined  in  phase  I 
of  the  project  is  a  nonlinear  continuum  theory  which  is  applicable  to  composite  lam¬ 
inates  that  are  composed  of  two  homogeneous  and  isotropic  constituents  within  the 
context  of  purely  mechanical  theory.  The  generalization  of  the  theory  to  include  the 
effect  of  temperature,  anisotropy  and  in  particular  stacking  sequence/fiber  orienta¬ 
tion  in  fiber  reinforced  composite  laminates  will  be  carried  out  in  this  task. 
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3.  Application  of  the  theory  to  various  practical  problems.  The  micro-macro  contin¬ 
uum  structure  built  in  the  kinematics  of  Cosserat  composite  theory  provides  an  ideal 
model  for  the  analysis  of  thick  composite  laminates  composed  of  a  large  number  of 
thin  plies.  Due  to  the  importance  of  the  interlaminar  stresses  and  edge  effects  at  the 
boundaries  of  composite  plates  and  shells  the  Cosserat  composite  theory  will  be 
utilized  for  the  analysis  of  composite  laminates  under  in-plane  and  out-of-plane 
loading  conditions.  The  existing  theories  that  address  the  interlaminar  stresses  are 
not  manageable  when  the  number  of  layers  exceeds  seven  or  eight  (Pagano,  1989). 
We  commented  earlier  on  the  significance  of  the  effect  of  discontinuity  in  compo¬ 
site  laminates  and  the  urgent  need  for  solutions  to  such  problems.  In  this  regard, 
the  analysis  of  an  initially  flat  composite  laminate  containing  a  circular  hole  will  be 
pursued  and  the  effects  of  the  discontinuity  will  be  studied. 

4.  Variational  formulation  of  the  Cosserat  composite  theory  suitable  for  finite  element 
discretization  of  composite  laminates.  This  part  includes  the  following  activities; 

•  formulation  of  field  equations  and  constitutive  relations  in  a  form  suitable  for 
finite  element  implementation  and  derivation  of  weak  forms  of  these  equations. 

•  linearization  of  the  weak  form  in  terms  of  linearized  strain  measures  and  finite 
element  interpolation  of  the  representative  Cosserat  surface  and  the  director 
field  and  approximation  of  various  stress  components. 

•  numerical  implementation  of  the  theory  and  evaluation  of  performance  of  the 
proposed  finite  element  model  in  solving  various  practical  problems. 

Our  efforts  during  the  year  1991  were  directed  toward  the  achievement  of  these  objectives. 


The  results  are  summarized  next. 
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U  Present  Status  of  the  Project 

The  results  of  our  research  efforts  during  the  last  year  are  presented  in  section  2.0 
through  section  7.0  of  this  report.  A  summary  of  the  contents  of  these  sections  is  presented  in 
the  following. 

In  section  2.0  the  kinematics  of  the  micro-  and  macro-structures  were  examined  and  the 
relationship  between  strain  measures  at  micro-  and  macro-levels  were  derived.  The  field 
equations  for  compos  .e  laminates  were  derived  through  a  direct  integration  of  field  equations 
of  classical  continuum  mechanics.  The  linearized  kinematic  measures  were  derived  in  the 
context  of  infinitesimal  deformation  and  the  relation  of  linear  strain  measures  with  displace¬ 
ment  vector  and  director  displacement  vector  were  obtained.  The  equations  of  motion  in  the 
linear  theory  were  derived  and  were  presented  for  both  curved  and  flat  geometries. 

Section  3.0  showed  the  derivation  of  constitutive  relations  for  composite  laminates.  A 
procedure  for  deriving  the  relation  between  composite  quantities  (i.e.,  composite  stress  tensor 
and  composite  couple  stress  tensor)  and  strain  measures  at  macro-level  were  presented.  The 
derivation  was  performed  for  a  bi-constituent  composite  laminate  and  the  constitutive  relations 
were  expressed  in  terms  of  material  constants  associated  with  every  individual  layer. 

Section  4.0  presented  the  complete  theory  for  linear  elastic  composite  laminates.  The 
relationship  between  the  displacement  vector  and  the  director  displacement  vector  was  derived 
based  on  the  geometrical  continuity  at  interfaces.  The  field  equations  were  derived  in  terms 
of  displacement  vector  and  it  was  shown  that  classical  continuum  theory  can  be  derived  from 
Cosserat  composite  theory  for  the  case  of  a  single  constituent.  The  theory  was  further 
simplified  for  bi-laminate  micro- structure  composed  of  isotropic  constituents.  Finally  the  con¬ 
stitutive  relations  for  composite  stress  tensor,  composite  couple  stress  tensor  and  interlaminar 
stress  vector  were  derived  in  terms  of  the  displacement  vector,  its  gradients  and  material  con¬ 
stants  of  individual  constituents. 
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Section  5.0  was  the  extension  of  the  theory  for  multi-constituent  composites.  The  micro- 
structure  or  representative  element  was  assumed  to  be  composed  of  several  constituents  which 
repeated  themselves  in  the  layering  direction.  The  development  of  this  section  is  particularly 
suited  for  fiber  reinforced  composites  where  the  fiber  direction  changes  in  the  stacking 
sequence  of  the  plies.  The  theory  was  simplified  for  the  case  of  isotropic  constituents. 

Section  6.0  presented  the  extension  of  the  theory  from  a  purely  mechanical  theory  to  a 
thermo-mechanical  theory.  In  this  section  composite  field  quantities  crresponding  to  the  heat 
flux  vector,  the  heat  supply  and  the  specific  entropy  of  classical  thermo-mechanical  theory 
were  introduced  and  the  equation  of  local  balance  of  energy  and  the  Clausius-Duhem  inequal¬ 
ity  were  derived  in  terms  of  these  composite  field  quantities. 

Section  7.0  presented  the  constitutive  relations  of  linear  thermo-eiasticity  for  composite 
laminates.  These  constitutive  relations  were  derived  for  the  composite  stress  tensor,  compo¬ 
site  couple  stress  tensor,  entropy,  heat  flux  vector  and  heat  flux  couple  vector.  The  develop¬ 
ments  of  this  section  were  parallel  to  those  of  section  4.0  and  a  set  of  coupled  thermo- 
mechanical  field  equations  in  terms  of  the  displacement  vector  and  the  temperature  were 
presented.  The  continuation  of  these  efforts  is  briefly  discussed  in  the  Future  Work  section. 

We  would  like  to  emphasize  that  because  of  the  very  complex  nature  of  composite 
materials,  any  general  continuum  theory  capable  of  representing  various  behavioral  aspects  of 
these  materials  will  be  complicated.  During  the  course  of  these  developments,  two  lines  of 
thought  regarding  the  interpretation  of  various  field  quantities  of  the  theory  were  considered. 
Accordingly,  some  of  the  earlier  developments  of  the  phase  I  were  re-examined  and  a  few 
modifications  regarding  the  definition  of  the  field  quantities  were  introduced.  This  process 
required  a  concentrated  effort  by  all  team  members,  but  it  had  its  own  rewards.  The  modeling 
and  the  mathematical  issues  were  clarified  and  a  better  understanding  of  various  features  of 
the  theory  was  achieved. 
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The  field  equations  for  composite  laminates,  as  presented  in  section  2.0  and  section  6.0, 
were  derived  directly  from  corresponding  field  equations  of  classical  continuum  mechanics. 
In  this  derivation  the  existence  of  a  micro- structure  in  the  form  of  a  Cosserat  surface  was 
assumed  for  the  composite  laminates.  A  parallel  development  of  the  theory  is  presented  in 
Volume  II  of  the  report  in  which  the  conservation  laws  of  a  single  micro-structure  (i.e.,  a 
Cosserat  surface)  were  integrated  across  the  layering  direction  of  the  composite  laminate. 
This  integration  process  provided  the  global  balance  laws  for  composite  laminates  and  the 
governing  field  equations  were  derived  from  these  global  balance  laws. 

Some  of  these  developments  were  presented  in  the  AFOSR  contractors  meeting  on 
Mechanics  of  Materials  in  October  1991  in  Dayton,  Ohio.  A  copy  of  this  presentation  is 
included  in  Volume  II.  Also  the  draft  of  a  technical  paper  presenting  the  field  equations  of 
composite  laminates  as  derived  from  conservation  laws  was  prepared  for  publication  in  the 
International  Journal  of  Engineering  Science.  A  copy  of  this  draft  report  is  also  included  in 
Volume  n  of  the  present  report.  The  research  efforts  presented  here  were  carried  out  by  Dr. 
M.  Panahandeh,  Dr.  G.  R.  Ghanimati,  Dr.  V.  Schricker  and  Dr.  M.  Mahzoon,  a  visiting  scho¬ 
lar  with  the  University  of  California,  Berkeley. 

1.4  Future  Work 

We  are  very  encouraged  with  the  developments  of  this  project.  These  developments  pro¬ 
vide  a  framework  for  the  modeling  and  analysis  of  various  composite  laminates  and  can  be 
extended  to  study  a  wide  range  of  technical  issues  associated  with  the  application  of  compo¬ 
site  laminates  such  as  nonlinear  material  behvior,  stability  analysis,  damage  growth  and  failure 
mechanisms,  to  name  a  few. 

The  continuation  of  the  present  work  to  achieve  the  technical  objectives  of  the  project,  as 
outlined  in  Section  1.2,  is  planned  for  this  year.  In  particular,  the  analysis  of  composite  lam¬ 
inates  under  in-plane  and  out-of-plane  loadings  and  the  numerical  solution  of  the  system  of 
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equations  derived  in  Section  4.0  is  given  immediate  priority.  Theoretical  development  of  the 
theory  for  cylindrical  and  spherical  geometries  is  underway.  The  derivation  of  constitutive 
relations  for  orthotropic  materials  presenting  fiber  reinforced  plies  with  various  fiber  orienta¬ 
tions  and  different  stacking  sequences  is  also  a  part  of  this  year’s  effort.  The  finite  element 
formulation  of  the  theory  and  development  of  a  composite  shell  element  will  be  followed 
immediately  after  the  completion  of  the  theoretical  developments. 
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2.0  MICRO-MACRO  CONTINUUM  MODEL  OF  COMPOSITE  LAMINATES 
2.1  Kinematics  of  Micro-  and  Macro-Structures 

Let  the  points  of  a  region  3^in  a  three  dimensional  Euclidean  space  be  referred  to  a  fixed 
right-handed  rectangular  Canesian  coordinate  system  x*  (i  =  1,2,3)  and  let  0‘  (i  =  1,2,3)  be  a  gen¬ 
eral  curvilinear  coordinate  system  defined  by  the  transformation  x*  =  x*(0j).  We  assume 

this  transformation  is  nonsingular  in  %,  Furthermore,  let  4  represent  the  coordinate  of  a  micro¬ 
structure  in  the  layering  direction  with  ^  =  0  corresponding  to  the  bottom  surface  of  the  micro- 
structure.  We  recall  that  a  convected  coordinate  system  is  normally  defined  in  relation  to  a  con¬ 
tinuous  body  and  moves  continuously  with  the  body  throughout  the  motion  of  the  body  from  one 
configuration  to  another. 

Throughout  this  work,  all  Latin  indices  (subscripts  or  superscripts)  take  the  values  1,2,3;  all 
Greek  indices  (subscripts  or  superscripts)  take  the  values  1,2  and  the  usual  summation  conven¬ 
tion  is  employed.  We  will  use  a  comma  for  partial  differentiation  with  respect  to  coordinates  0“ 
and  a  superposed  dot  for  material  time  derivative,  i.e.,  differentiation  with  respect  to  time  hold¬ 
ing  the  material  coordinates  fixed.  Also,  we  use  a  vertical  bar  ( I )  for  covariant  differentiation. 
In  what  follows,  when  there  is  a  possibility  of  confusion,  quantities  which  represent  the  same 
physical/gcometrical  concepts  will  be  denoted  by  the  same  symbol  but  with  an  added  asterisk  (*) 
for  classical  three  dimensional  continuum  mechanics  and  no  addition  for  composite  laminate 
(macro-structure).  For  example,  the  mass  densities  of  a  body  in  the  contexts  of  the  classical  con¬ 
tinuum  mechanics,  and  the  composite  laminate  (macro-structure)  will  be  denoted  by  p*  and  p, 
respectively. 

The  micro-macro  continuum  model  of  a  composite  laminate  is  illustrated  in  Figures  1  and 
2.  Figure  1  shows  a  typical  composite  laminate  (only  three  micro-structures  are  shown  in  this 
figure).  Figure  2  shows  a  shell-like  micro- structure  with  its  associated  coordinates.  This  micro¬ 
structure  is  composed  of  two  constitutents  and  can  be  generalized  for  cases  of  multi-constituent 
composites. 
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micro  structure 


Figure  1 

A  composite  laminate  consisting  of  alternating 
layers  of  two  materials 


Figure  2 

A  SHELL-LIKE  MICRO-STRUCTURE  (REPRESENTATIVE  ELEMENT) 
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We  begin  the  development  of  the  kinemadcal  results  by  assuming  that  the  position  vector 
of  a  particle  P*  of  a  representative  element  (k*  micro-structure),  i.e.,  p*(0“,0^\4,t)  in  the 
present  configuration  has  the  form 

p*  =  r(0“.03«,t)  +  4d(0“.0^^\t)  (k  =  (2.1) 

where  r  is  the  position  vector  for  the  surface  4  =  0  **  is  the  director  field.  0^^\  at  this 

point,  is  an  identifier  for  the  k*  micro-structure.  Greek  super-  or  subscripts  will  assume 
values  of  1  and  2  only.  The  dual  of  (2.1)  in  a  reference  configuration  is  given  by 

P*  =  R*(0“,03<*'))  -I-  4D(0“,03Oc>)  (2.2) 

If  the  reference  configuration  is  taken  to  be  the  initial  configuration  at  time  t  =  0,  we  obtain 

p*(0o,03(k)  ^,0)  =  r(ea^e30'\O)  4d(0“  0^^>,O) 

=  R(0“,03<''>)  -I-  4D(0“  0^*))  =  P*(0“,03<^>,4)  (2.3) 

The  velocity  vector  v*  of  the  three-dimensional  shell-like  micro-structure  at  time  t  is 
given  by 


* 

V 


=  p*(0“,03«,4.t) 


(2.4) 


where  a  superposed  dot  denotes  the  material  time  derivative,  holding  0°  and  4  fixed.  Rom 
(2.1)  and  (2.4)  we  obtain 


where 


v*  =  v  +  4w 


(2.5) 
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V  =  r 


w  =  d 


(2.6) 


The  base  vectors  for  the  micro-  and  macro-structures  are  denoted  by  g*  and  gj,  respec¬ 
tively,  and  we  have 


^p*  9p* 

Using  (2.1)  and  (2.7)  we  obtain  the  following  relations  between  g*  and  gj 


(2.7) 


=  ga  +  a 
83  =  83  =  «• 

where  (  )„  denotes  partial  differentiation  with  respect  to  6". 


(2.8) 


By  a  smoothing  assumption  we  suggest  the  existence  of  continuous  vector  functions 
gi(0“,0^)  for  the  macro-structure  with  the  following  property 


8j(0®  ~  gi(0“.e3<^>)  (2.9) 

where  g(0®,0^^^)  are  defined  according  to  (2.7)2-  ^  similar  smoothing  assumption  is  also 
made  for  the  director  d  which  we  like  to  attach  to  every  point  of  the  macro-structure.  Based 
on  the  smoothing  assumptions  we  can  write  (2.8)]  as  follows 


8a  ~  8a  ^8k  (3  a) 

where  (  )  stands  for  the  Christoffel  symbol  of  the  second  kind  and  is  defined  as 


(2.10) 


(j  M  =  g«[3aoJ«  =  i 


^83a- 

80j  ^ 
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The  following  relations  can  also  be  derived  between  the  components  of  metric  tensors 
gij  =  g*  ■  g*  and  gij  =  gi  •  gj 


fep  *  gap  +  ^1(3  ^  a}gpk  +  {3  ''  p)gakl  +  a)  I3  ^  p)gkj 

g^  =  go3  +  ^3 a)gk3  (2.11) 

S33  =  g33 

which  after  simplification  and  linearization  in  terms  of  ^  reduce  to 


gap  =  gap  +  4gap,3 

ga3  =  ga3  +  Y  4g33.a  (2.12) 


g33  -  g33 

The  determinants  of  metric  tensors  gy  and  gy  arc  also  related  according  to  the  following 
relation 


g*  =  g  +  ^ 


(2.13) 


where 


g*  =  det(gy)  ,  g  =  det(gy) 


gll.3glZ3g33,l 

gll  gl2  gl3 

A  = 

gl2  g22  g23 

+ 

gl23g2Z3g33a 

gl3  g23  g33 

gl3  g23  g33 

and  the  final  result  has  been  linearized  in  terms  of 

We  recall  that  the  director  d  is  the  same  as  g3  and  therefore  when  referred  to  the  base 
vectors  gj  it  has  only  one  non-zero  component,  namely  d^  =  1,  so  we  can  write 
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d  =  d‘gi  ,  d“  =  0  ,  <P  =  \ 
di^gijCjj  ,  di  =  gi3  (i=l,2.3) 


(2.15) 


where  dj  and  d‘  denote  the  covariant  and  contravariant  components  of  d  referred  to  g‘  and  gj, 
respectively.  The  gradient  of  the  director  d  may  be  obtained  as  follows 


di  =  g3.i  =  {3  iJgk  =  d^igk 

The  vertical  bar  (  I  )  denotes  covariant  differentiation  with  respect  to  gjj.  For  convenience  we 
introduce  the  notations 


=  gi  •  =  diij 

X)  =  t  •  d  j  =  dSj 

From  (2.17)  it  is  clear  that 

=  g%j 

Making  use  of  (2.17)  and  (2.16)  we  have: 

^■j=dSj={3^) 

=  gki^'^j  =  I3j4] 

Consider  now  the  velocity  vector  v  which  can  be  written  in  the  form 


(2.17) 


(2.18) 


(2.19) 


v  =  v*gi  =  Vjg'  (2.20) 

Again  we  make  a  smoothing  assumption  for  the  existence  of  the  vector  function  v(0“,6^)  such 
that  v(0“,0^)ie3^3(k)  =  r(0“,0^^^)  after  which  we  can  define  the  gradient  of  the  velocity  field 
and  we  have 


v.i  =  (V^gj).i  =  vj|igj 


(2.21) 
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We  now  introduce  the  notations 


From  (2.22)  it  is  clear  that 


vy  =  gi  •  Vj  =  vjij 
v‘j  =  g‘  •  Vj-  =  vSj 


(2.22) 


v‘j  =  g^'kj 

V.i  =  Vygj  =  vjjgj 


(2.23) 


We  observe  that  both  Xy  and  vy  represent  the  covariant  derivative  of  vector  components  anH 
hence  transform  as  components  of  second  order  covariant  tensors. 


We  may  decompose  vy  into  its  symmetric  and  skew-symmetric  pans,  i.e.. 


Vy  =  +  Vjyj  *  T\ij  +  CO-y 

^ij  =  ''(ij)  =  7  (VV  =  (2.24) 

=  '"[ijl  =  7  (''ir''ji)  =  -  «ji 

Also  in  view  of  (2.6),  (2.7)2  (2.24)i  we  may  express  gj  in  the  form 


ga  =  V.a  =  (Tlka  + 
g3  =  d  S  w  =  Wfcg*"  =  W^g,^ 

The  gradient  of  the  director  velocity  in  0“  direction  is  obtained  by  writing 


(2.25) 


W.O  =  <*.a  =  (d.„)  =:  X^gk 


=  ^  gk  +  ^gk  =  ^gk  +  ^gp  + 


=  ^gk  +  ^|(Tlkp-«Oicp)g‘‘  +  ^aW‘‘gk 

=  [^  +  >^(ti||-kdJ)  +  ^W^]gk 


(2.26) 

The  dual  of  the  above  expressions  in  the  reference  configuration  can  be  written  easily  by 
substituting  appropriate  capital  leners  for  small  letters. 

We  now  introduce  relative  kinematic  measures  y-j  and  Ky  such  that 

■Kj  =  I  (*ij  -  G|j)  =  Jf  (2-27) 

Kjj  =  Xij  -  Aij  (2J8) 

where 

Gy  *  Gi  •  Gj  (2.29) 

Aij  =  t3j.i]  =  -i  (Gay  +  Gji,3-G3i4)  (2.30) 

Making  use  of  (2.12)  and  similar  expressions  for  the  reference  configuration  we  can 
relate  relative  kinematic  measures  of  the  micro-structure  as  follows 

Vi  =  Ypo  =  (gi  -  G^)  =  Y  [(gap  +  ^  gap^)  -  (Gap  +  ^  Gopj)] 

=  Yop  +  ^  ^(gaP3  ”  ^op^) 

~  Yop  ~2  ^^^P  ^p<x)  (2.31) 

Y^  -  Yto  =  "J  I(go3  *2  ^  ®33.a)  “  (Ga3  ^  ^  G33  a)] 

=  Ya3  +  4  K3a 


(2.32) 
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Y33  =  %3 

In  obtaining  the  above  results  we  have  noted  that 


X(jtp  =  [3p,a]  -  —  (g3tt^3  -  g3p  „)  +  -j  3 
^3a  =  [3cx,3]  =  —  §33  „ 

Kap  +  Kpoj  =  gjjjpj  -  Gojp  3 

and  we  have  linearized  the  result  in  terms  of  4. 


(2.33) 


(2.34) 
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2.2  Basic  Fidd  Equations  for  Micro*  and  Macro*Structures 


The  three-dimensional  equations  of  motion  in  classical  continuum  mechanics  are 
recorded  here  for  the  k***  representative  element  (micro-structure)  in  the  present  configuration 

=  0  (2.35) 

T*\i  +  p*bV^^  =  p*v*g**^  (2.36) 

g*xT*‘  =  0  (2.37) 

where 

t*  =  g*-i/2T*in  •  ^  x*'  =  g**V«^*  (2.38) 

The  argument  of  all  starred  functions  recorded  above  is  (0“,0^^\^,t)  and  the  equations  are 

written  for  each  and  every  representative  element  (k  s  i,2 . n)  which  is  assumed  to  repeat 

itself  in  the  present  model. 

Now  introduce  the  following  quantities  for  each  micro-structure: 

Composite  Stress  Vector  T‘: 


S2  0 


Composite  Stress  Couple  Vector  S“: 


sa(ea  03(k)  t)|:±  f  ^T*®(0“,03O'>4.t)d^ 
S2  0 


Composite  Mass  Density  p: 


(2.39) 


(2.40) 


(2.41) 
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,  ^ 

pgi^(z“)i:-—  f  (^)‘*pV*'^d4  (as:l,2) 

^2  0 


(2.42) 


Composite  Body  Force  Density  b: 


J  pVg*^^d4 


(2.43) 


Composite  Body  Couple  Density  c: 


pg>“ci-i/pVg*>«5d4 


(2.44) 


The  quantities  on  the  left-hand  side  of  equations  (2.39)-(2,44)  are  discrete  in  terms  of  the  vari¬ 
able  6^^^  which  are  made  continuous  by  smoothing  assumptions.  The  composite  mass  density 
Po  in  the  reference  configuration  is  also  defined  as  follows: 


PoC"’  =  4-  j  PoV'^dS 
^2  0 


(2.45) 


where  pj  is  the  mass  density  of  the  micro-strucnire  in  the  reference  configuration.  Since 
p*g*i/2  _  p*(3*i/2^  jjjg  continuity  equation  for  the  macro-structure  is  readily  seen  to  be 


(2.46) 


Now  consider  equation  (2.36)  and  first  divide  it  by  ^2  and  then  integrate  with  respect  to 
4  from  0  to  ^2  obtain  the  equation  for  balance  of  linear  momentum  for  the  macro-structure 


=  ~  J  p*(v  +  4w)g*''^4 


(2.47) 


Each  term  in  the  above  equation  can  be  represented  in  terms  of  the  quantities  defined  in 
(2.39)-(2.44)  except  the  second  term  which  is  the  difference  between  interlaminar  stresses 
above  and  below  the  representative  element  divided  by  its  thickness  ^2 

T-  J  =  -r  -  T*3(0“,e3O'>,t)]  (2.43) 

S2  0  ^2 

Now  we  postulate  the  existence  of  the  continuous  vector  function  o(0“,0^,t)  whose  values  at 
03  _  03(k)  are  the  same  as  interlaminar  stresses  T*^(0“,0^dc)  and  further  approximate  (2.43) 
as  the  gradient  of  this  function  in  the  0^  direction.  With  this  in  mind  we  write  (2.47)  as 

T“a  +  =  pg*'^(v  +  z^w)  (2.44) 

O0'^ 

To  obtain  the  equation  for  balance  of  director  momentum,  (2.36)  is  multiplied  by 
integrated  from  0  to  ^2  divided  by  ^2  to  got 

, 

=  T”  f  P*g*^'^(^v  +  ^^)d^  (2.50) 

S2  0 

Again  the  second  term  in  the  above  equation  can  be  wrinen  as 

^  3rr*3  ^ 

5  d?  =  i  J  T'’d5  =  o  -  T’  (2.51) 

As  a  result  we  have 

S“o  +  o  -  +  pg*^c  =  pg’^(z*v  +  z^)  (2.52) 


which  is  the  equation  for  balance  of  director  momentum. 
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Ncxt,  wc  consider  (2.37),  divide  it  by  ^2  integrate  with  respect  to  ^  ft®™  ^=0  to 
^=^2  2nd  making  use  of  (2.8)i^  we  get 


■T  I  T*“  +  g3  X  T*3)d^  =  0 

^2  0 

1  1  ^ 

J  (ga  +  ^da)>^T*“d^  +  -^  JdxT*3d^  =  0 


(2.53) 


and  substituting  from  (2.39)  and  (2.40)  we  obtain 


g„  X  T“  +  d  „  X  S“  +  d  X  =  0 


which  can  also  be  written  as 


(2.54) 


Si  X  T‘  +  da  X  S“  =  0 


(2.55) 


This  is  the  balance  of  moment  of  momentum  for  the  macro-structure. 


Now  we  proceed  to  obtain  an  expression  for  the  specific  mechanical  energy.  Such  an 
expression  for  each  micro-structure  can  be  written  as 


'T>*i  • 

p  g  =  T  ‘  •  V  i 


(2.56) 


First,  using  (2.5)  we  write  this  equation  as 


p*g*i/2^*  _  j*a  .  ^  (V  +  ^w) 


=  T*“  •  v,„  +  ^T*“  •  w.a  +  • 


(2.57) 


Dividing  (2.57)  by  ^2  and  integrating  with  respect  to  ^  from  ^  =  0  to  ^  =  ^2  will  result  in 


j-  j  P'g‘‘'^’d5  =  /  T'«d4  ■  V  a  +  -^  /  5T‘“d5  ■  "  o  i  ■ 


(2.58) 


BASE 


We  now  define  composite  specific  mechanical  energy  for  the  representative  element  by 


pg*^e  =  T-  f  (2.59) 

^2  0 

From  this  definition,  the  equation  of  continuity  and  other  definitions  (2.39)  through  (2.44)  for 
composite  quantities,  (2.58)  can  be  rewritten  as 

pgi/2£  =  .  V  „  +  S“  •  w  o  +  '  w  (2.60) 

Since  v  =  r,  v  ^  =  (r)  „  =  ^  =  go  and  w  =  d  =  fo,  we  can  further  reduce  (2.60)  to 


pgi/2g  =  X'  .  g.  +  s“  •  w 


.a 


(2.61) 


which  is  the  appropriate  expression  for  the  specific  mechanical  energy  of  the  macro-structure. 


T“a  +  +  pg^^b  =  pg^'2(v  +  z^w) 

(2.62) 

S“  <ji  +  o  -  +  pg'^c  =  pg^^(z*v  +  z^w) 

(2.63) 

gj  X  T’  +  X  S“  =  0 

(2.64) 

And  also  the  following  expression  was  derived  for  the  specific  mechanical  energy 


pgi/2e  =  T‘-gi  +  S“  - w,a  (2.65) 

By  referring  various  vector  quantities  to  the  base  gj  we  would  like  to  write  the  above 


equations  in  component  form.  First  write 

Ti  =  gi/2  xijg.  (2.66) 

a  =  o'gj  (2.67) 

S“  =  g»^  S“jgj  (2.68) 

b  =  b*gj  ,  c  =  (jjgj  (2.69) 


where  and  8“^  are  contravariant  components  of  composite  stress  tensor  and  composite  cou¬ 
ple  stress  tensor,  respectively,  o'  is  the  interlaminar  stress.  Now  substitute  in  (2.62)  and 
obtain 


(gl/2  ^  ^  ^  pgl'25jg.  =  pgl/2(^  + 

'C“^).agj  +  g^'^  t“Mj  o)gk  +  Ij  S)gk  +  Pg''^b>gj 


or 


(gl/2  +  gl/2  ^ak{^  j  al  +  (k  ^  3}  +  P 

=  P  g^^  (^  +  Z^'^) 


(2.70) 


Equation  (2.63)  reduces  to 


or 


S®Jgp  „  +  o>g.  _  gl^  X^jg.  +  pg»'2cjg. 

=  pg''^(z^^  +  Z^w’)gj 


(gl/2  s“j)  ^  +  g^^  S®*  {,£  j  +  d  -  g*^  T^j  +  pg^'^cj  =  p  g^'^(z'v^  +  z^w*)  (2.71) 

Equation  (2.64)  can  be  rewritten  as 


gi  X  (gJ/^  xUgj)  +  X'ogi  X  S®jgj)  =  0 

or 

gi/2(^ij  +  Xi„S®j)gi  X  gj  =  0  (2.72) 

since  g  ^  0,  gj  x  gj  =  e^j^g''  and  Ejjij  is  skew-symmetric  we  conclude  that  the  quantity  in 
parentheses  in  (2.72)  must  be  symmetric  in  i  and  j.  As  a  result,  the  conser/ation  of  angular 
momentum  in  component  form  is  the  symmetry  of  defined  by 


•pj 

(2.73) 

•pij  _  'pi 

(2.74) 

The  expression  (2.65)  for  the  specific  mechanical  energy  can  also  be  written  as 
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pg^  *  t^gj  •  gi  +  S“jgj '  w'logi 

-  gl/2(t®ig.  .  ^  +  S“jwj|„  +  T^jgj  •  fo) 

-  gi^(t“jgj  •  v„  +  T^jgj  •  w  +  S®^,„) 

=  g^'^(t“jvj,Q  +  t^jwj  +  S“jwj,a) 

We  have  now  the  component  form  of  the  expression  for  mechanical  power 

P  i  pe  =  t^Vj ,  „  +  t^Jwj  +  S“jwj  I  a  (2.75) 

An  alternative  form  for  mechanical  energy  expression  is  derived  in  which  the  rates  of 
relative  Icinematic  measures  will  appear.  Using  (2.25)]  and  (2.26),  we  rewrite  (2.75)  as 

P  =  pe  =  t“j(Tij„  +  (Ojo)  +  t^jwj  +  SPj[Xjp  +  (Tija  +  a)j„)  +  X|wj] 

=  (t“j  +  S%“)ilj„  +  (x“j  +  sWXp“)(Dje,  +  sW^p 

+  (T^j  +  X|sPj)Wj  (2.76) 

Recalling  (2.73)  and  using  symmetry  of  and  skew-symmetry  of  Wy  we  can  write  (2.76)  as 

P  =  T«Mja  +  sPJXjp  +  T^iwj  +  T«3co3„  (2.77) 


By  (2.25)i  we  have 


go  •  gp  =  ^pa  +  <«>pa  >  gp  •  go  =  "nop  +  0)ap 


(2.78) 


Therefore 


gop  =  go  '  gp  go  '  gp  =  2t1oP 


(2.79) 


1 

^op  “  2  ~ 


(2.80) 
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In  the  last  result  we  have  used  the  definition  of  Yap  from  (2.27).  By  (2.25)  12  we  have 


TherefOTe 


ga  ■  83  -  %o  +  0)30 
83  •  ga  =  Wa 


(2.81) 

(2.82) 


ga3  =  ga  •  83  +  83  •  ga  =  ^30  +  tOjct  +  W„  (2.83) 

'n3o  ~  goG  ~  (0)3o[  +  Wjjj)  =  2Ya3  ■”  (-^^a  (2.84) 

Again  by  (2.25)2  and  (2.27) 


=  Kjp  (2.87) 

Substituting  from  (2.80),  (2.84),  (2.86)  and  (2.87)  in  (2.77)  we  get 

p  =  T^oS  +  T«5(2r.3  -  ««3.  +  Wa))  +  s"^!)  +  T’«w„  + (2.88) 
which  is  simplified  to 

P  =  T“^oP  +  2T^Ya3  +  T^S3  +  sPJkjp  (2.89) 

If  symmetry  of  and  y-j  is  considered  we  can  further  simplify  (2.89)  and  obtain 

P  =  +  SPjii^p  (2.90) 

or 
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2.4  General  Constitutive  Assumption  for  Elastic  Composite 

At  this  point  we  postulate  the  existence  of  specific  internal  energy  in  purely  mechanical 
theory  which  depends  on  relative  kinematic  measures  yy  and  as  defined  in  (2.27)  and 
(2.28) 


¥  =  ¥(Yy.Kja) 

P  =  PV 

By  usual  procedures  we  obtain  from  (2.91),  (2.92)  and  (2.93) 


(2.92) 

(2.93) 


S“j  =  p  ^  (2.95) 

dKja 

Now  the  composite  stress  vector  T'  and  the  composite  couple  stress  vector  S“  from  (2.66)  and 
(2.68)  will  be 


S«  =  pgi/2(^)g  (2.97) 

The  coefficient  pg*^  can  be  replaced  by  poG*^  by  taking  advantage  of  the  continuity  equa¬ 
tion.  Note  that  by  these  constitutive  relations  for  P  and  S°  the  balance  of  moment  of 
nromentum  is  identically  satisfied. 
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15  Linearized  Kinematics 

For  linearized  kinematics  let 


r(0a  e3(k)  t)  =  R(e«  e^<«)  +  eu(e“,e3*\t) 

(2.98) 

d(e“,0^^>,t)  =  D(e“,0^<^^  +  e5(0“,0^(*'\t) 

(2.99) 

•  • 

v  =  r  =  eu  ,w  =  d  =  e5 

(2.100) 

where  e  is  a  non-dimensional  parameter.  The  motion  of  the  macro-structure  describes 
infinitesimal  deformation  if  the  magnitude  of  the  gradient  of  the  displacement  vector  eu  and 
the  magnitude  of  the  director  displacement  vector  e5  are  of  the  order  of  e  «  1  such  that  in 
the  following  developments  we  can  only  retain  terms  which  are  linear  in  e.  The  base  vectors 
gi  are  found  from  (2.7)2 


ga  =  R  a  +  eu  a  (2.101) 

g3  =  d  =  D  +  e5  (2.102) 

The  corresponding  vectors  in  reference  configuration  are; 

Ga  =  Ra  .  G3  =  D  (2.103) 

We  now  proceed  to  obtain  the  relative  kinematic  measures  Yij  ^d  KjQ.  Using  (2.103)2 
(2.101)  together  with  the  definition  of  g^^  and  Gap  we  write 

Sop  ~  (Gq  +  eu  a)  •  (Gp  +  eu  p)  =  Gap  +  e(G„  •  u  p  +  u  „  •  Gp)  +  O(e^)  (2.104) 

where  O(e^)  denotes  terms  of  order  e^  in  displacement  gradient,  where 

Gq  ‘  “p  “a  ■  Gp  =  Gq  •  iPipgj  +  uj|„gj  •  Gp 

*  Ga  •  u7,pgy  +  Ga  •  U^ipg3  +  u''^|aBy  '  Gp  +  u^|ag3  •  Gp 


BASE 


-31  - 


=  u^ipGot  •  (G^  +  eu.y)  +  u^ipGa  •  (D  +  e5)  +  u^|<x(Gy  +  £11^)  ■  Gp 

+  u^(g(D  +  eS)  '  Gp  (2.105) 

Retaining  terms  which  are  of  the  order  of  unity  in  (2.105)  and  substituting  the  result  in 
(2.104)  we  find 

Yap  =  -J  (gop  -  Gap)  =  J  (Uolp  +  Upla)  +  J  (ufpDo  +  U?aDp)  (2.106) 

Here  covariant  differentiation  is  supposed  to  be  performed  with  respect  to  the  metric  Gjj  of 
the  reference  configuration  and  instead  of  £u  we  have  used  u  with  the  same  assumptions  made 
for  linearization.  Similarly  we  can  write 

go3  =  (Gjt  +  EU  <x)  •  (G3  +  £6)  =  Ga3  +  E(G(jt  •  5  +  U  ^  •  G3)  +  0(£^) 

=  Ga3  +  £(5„  +  U3,a)  +  0(e^)  (2.107) 

Again  using  5  instead  of  eS  with  the  same  interpretation  we  obtain 

Ya3  ~  Y3a  ”  (Sa3  “  G03)  =  ~  (Sq  +  U3|5[)  (2.108) 

To  find  Y33  we  write 

833  ~  (G3  +  £6)  ’  (G3  +  e5)  =  G33  +  2E53  +  O(e^)  (2.109) 

733  =  "J  (g33  "■  G33)  =  83  (2.110) 

As  for  the  measures  Kja  we  proceed  as  follows 

^  =  ga  •  <*.p  =  (Go+eUo)  •  (D,p^■E5p)  =  Aotp  +  e(uj,agj  •  Dp  +  G^  •  5j,pgj)  +  0(e^)  (2.111) 
where: 
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*  Dp  *  uT^ioCGy  +  eu,Y)  •  Dp  +  u^iaCGj  +  e6)  •  Dp 

=  ■•■  'i^ia^ap  ■*■  G(e) 

Go  '  ^ip8j  =  ^ipGo  *  (Gy  +  eu y)  +  S^ipGjx  •  (G3  +  e8) 

=  5a,p  +  83,pDa  +  0(e) 

Substituting  these  results  in  (2.101)  and  using  the  definition  of  Kqp  we  get 

Kqp  =  Xqp  -  Aqp  =  uj|aAjp  +  60, p  +  5^ipDo 

Now  we  obtain  an  expression  for  K3Q 

^30  =  83  •  <*.a  =  (G3  +  e5)  •  (Do  +  e5 o) 

^3a  =  A30  +  e(G3  •  5_a  +  5  •  D  „)  +  O(e^) 

We  simplify  each  term  separately 

G3  •  5.0  =  G3  •  (S^iogj)  =  5T^ioG3  •  (Gy  +  Eu,y) 

+  53|aG3  •  (D  +  e8)  =  S^ioDy  +  53,oD3  +  0(e) 
=  5j,oDj  +  0(e) 

S  •  D.O  =  (5igj)  •  (DSoGk)  =  (8igj)  •  d3,oG3 
=  Aa(5^gy  +  5^83)  •  G3 
=  A  3[5Y(Gy  +  eu,y)  +  8^(G3  +  e5)]  •  D 
=  0(e)  +  AjS^Dy  +  A  3530^  =  A^SJDj  +  0(e) 


(2.112) 


(2.113) 


(2.114) 


(2.115) 


However,  since  D“  =  0  and  =  1 
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5jDj  =  5jDj  =  63  (2.116) 

Substituting  from  (2.114),  (2.115)  and  (2.116)  in  (2.113)  and  using  previous  notation  for  5  wc 
obtain 

*^3a  ~  ^3a  “  ■^3a  ~  ^ta^j  “  ^lo  "*■  (2.117) 

To  recapitulate  the  relative  kinematic  measures  in  linear  theory  are: 

Yap  =  ^  (Uoip  +  Upio  U^laDp  +  U^ipDa) 

Ya3  =  Y3a  =  Y  (“31a  +  ^a) 

733  =  ^3  (2.118) 

Kop  =  ApUjia  +  5aip  +  S^ipDa 

•^Sa  “  ^la  AJ^ 

For  a  composite  with  initially  flat  plates  we  can  always  choose  our  base  vectors  such  that 
Gy  =  G“  =  5ij  and  as  a  result  =  Dq  =  0  and  if  we  confine  ourselves  to  small  deformations, 
then  all  Qiristoffel  symbols  vanish  and  covariant  differentiations  reduce  to  partial 
differentiations  and  equations  (2.1 18)  for  relative  kinematic  measures  will  funher  reduce  to 

Yap  =  y  (“a.p  +  “p.a) 

Ya3  =  y  (“3.0  +  5a) 

(2.119) 

Y33  =  53 

KqP  =  5o,p  ,  K3<x  =  63  a 

In  writing  these  relations  it  has  been  noted  that  =  0,  D3  =  1  and  Ai  s  0. 
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It  is  also  desirable  to  find  a  relation  between  g  and  G,  determinants  of  metric  tensors  in 
present  and  reference  configurations,  in  the  linear  theory.  We  recall  the  following  relations 

gl/2  =  gi  X  g2  •  g3  ,  =  Gi  X  G2  •  G3 

gi  X  gj  =  (Gj  +  eu,j)  X  (G2  +  eu_2) 

=  Gj  X  G2  +  e(ui  X  G2  +  Gj  X  u^)  +  O(e^) 

[818283!  =  [Gi  X  G2  +  e(ui  X  G2  +  G,  X  u^)]  •  IG3  +  eS]  +  O(e^) 

—  G^^  +  £[G]  X  G2  ‘  6  +  ii|  '  G2  X  G3  +  G3  X  G]  '  ii_2]  +  0(£^) 

=  G*^  +  £G’^[G3  •  5  +  u,i  •  G*  +  G2  •  U.2]  +  0(£2) 

Retaining  terms  which  arc  of  the  order  of  £  and  using  the  previous  notations  for  u  and  5  we 


(X)i/2  =  1  +  53  +  u«,„  (2.120) 

G 

Now  the  equation  for  balance  of  mass  will  readily  reduce  to 

Po  =  =  P(1  +  8^  +  ““i«)  «121) 

and  since  in  linear  theory  displacement  vector  u  and  director  displacement  5  satisfy  linearity 
assumptions  we  obtain 

P  ”  Po(i  ~  5^  “  u*ia)  (2.122) 
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2.6  Linearized  Field  Equations 

We  use  pertinent  results  from  linear  kinematics  and  usual  procedures  for  linearization  to 
write  the  held  equations  in  linear  theory.  It  should  be  recalled  in  such  analysis  that  g  is 
replaced  by  G,  p  by  and  Christoffel  symbols  are  calculated  with  respect  to  Gy.  By  omit¬ 
ting  the  details,  the  linear  version  of  the  equations  of  motion  are  recorded  here 

(G*'^  x“j)„  +  j  ai  +  ^  3}  +  p„  bi  =  Po  G^'^  (ip  +  (2.123) 

(G^'^  S“%  +  5“*^  (k  J  a)  +  <^  -  +  Po  cJ  =  =  Po  G''^  (z^iP  -I-  (2.124) 

T‘j  =  t'J  +  AiS“j  =  Tj‘  (2.125) 

For  a  composite  with  initially  flat  plies  further  simpliflcadon  can  be  made.  As  men¬ 
tioned  earlier,  G  =  1  and  all  Christoffel  symbols  vanish  identically.  The  resulting  balance 
equations  for  such  a  situation  will  be 


+  <^.3  +  Pob*  =  Po(u^  + 

(2.126) 

S“^,a  +  Cfj  -  +  PoCj  =  Po(z'iP  +  Z^^) 

(2.127) 

T'i  =  x'j  = 

(2.128) 
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3.0  CONSTITUTIVE  RELATIONS  FOR  LINEAR  ELASTICITY 


For  the  representative  micro-structure  let 


=  .  a  =1.2 

where 


(3.1) 


r 


0  <  ^  < 
4i<^<42 


(3.2) 


and  C(^^  (a  -  1,2)  are  material  constants  in  associated  layers.  Now  we  proceed  to  calculate 
T‘  and  S“  defined  in  (2.39)  and  (2.40).  First  we  recall  that  T*'  = 

g*i/2  _  gi/2  (2  ^  i  =  gy  +  g3*  =  g3  and  for  brevity  we  omit  the  index  a  in  rela- 

dons  (3.1)  and  (3.2) 


.  ^2  ^ 

Substitute  from  (2.31),  (2.32)  and  (2.33)  in  the  above  relations  and  get 


T  =  j  g'l'^  (c‘i»»Ya()  +  2c'i"Ya3  + 

¥  I  +  K3ac‘X^Kgj'<l4 

=  TW J  0*1“  g*'-^  g;  d4  +  {  (K.^  +  Kp,)  J  g-'^  gj- 5  d5 
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+  K3a  I  g*  ^  (3.3) 

We  calculate  each  integral  separately,  noting  that 

g”''^  h  =  g''^  tg,  +  4(-^  g,  +  <l,r)  +  (3.4) 

g*i«  g*  =  gW  (gj  +  A  |g3)  =  gin  (1  +  A  5)g3  (3.5) 

The  first  integral  in  (3.3)  is 

J  g-in  gi*  d?  =  A  J  (cHUg-in  g; g-in  g.)^^  (3  j) 

The  first  term  of  (3.6)  from  (3.4)  is  equal  to 

8r  1  I  ^  i  }  ^2  c‘TlWd^  (3.7) 

and  its  second  term  can  be  written  as,  from  (3.5), 

1  ^  A  1  ^ 

Combining  (3.7)  and  (3.8)  we  rewrite  (3.6)  as 


gin  g.  A  I  g.  j  gin  ^  j 

(3.9) 

The  second  integral  in  (3.3)  is 
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X  j  c»jo^  g*i'2  gj*  ^  J  g**'^  gy*  +  (3.10) 

^2  0  S2  0 

The  first  term  of  (3.10)  from  (3.4)  is 

1  ^  .  A  1  . 

gy  J  gy  +  d  y)  j  ^^c‘'>“Pd^ 

»•"> 

and  its  second  term  from  (3.5)  is 

Combining  (3.11)  and  (3.12)  we  rewrite  (3.10)  as 

^  ^2  ^ 

Si  f  /  + TT?  *i  f  /  *■''  ^  r  / 

S2  0  2g  S2  0  ^2  0 

The  third  integral  in  (3.3)  is 

1  j  =  g'-^  gj  -f  J  ^ciio’d?  + 

S2  0  ^2  0  ^8  ^2  0 

+  g>'2  d,  j  5Vl"d4  +  ^  ".T  -f  J  ^’<='^‘‘5  <3-“'> 

^2  0  ^8  ^2  0 

The  last  result  was  written  by  noting  the  development  in  (2.13).  The  results  in  (3.9),  (3.13) 
and  (3.14)  can  further  be  simplified  by  recalling  (2.9)3,  namely  dy=  and  using  the  fol¬ 
lowing  definitions  and  results. 
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Define 


m  =  <  1 


(3.15) 


a  = 


ai  0  <  ^  <  4i 

I  aa  <  4  <  42 


Then 


/5'‘ad^  =  A-t™,,  +  (^-m)ad 


k^-1 


m* 


Now  (3.9)  is  equal  to 


(3.16) 


(3.17) 


gj{mcp'  +  (l-m)cp^  +  ~ 

+  ~  [mep^  ■'■  ^  +  (-—  -  m)c^f^]] 

=  gj{(l  +  ^ 


.  1+na  4i(l-m)  , 

+  (l-m)(l  +  — - )cp'  +  — - - Xi(l  +  m 

4m  g  ‘‘  2m  ” 


.  4iA  1+m+m^^  iMf/, 


(3.18) 


(3.13)  can  also  be  written  as 


gjdmcpP  +  (-i-  -  mlcpP]  ^  +  A  .k  [mcij“P  +  (A 

m  2  2g  3  xir 


[mcl^^  +  (A-  -  m)c2^]) 

m-* 


g“  (j  Ky  +  ■^)inci>^  T  I 


+  (1  +  m  + -1^  + -^[4  - 

2m  3g  m  3 


m 


+  4  — (^  -  m)lcii»fl)  (3.19) 

o  g  nr 

The  expression  for  (3.14)  is  exactly  the  same  as  (3.19)  except  that  (oc,p)  in  (3.19)  should  be 
replaced  by  (a,3).  These  results  should  be  incorporated  in  (3.3)  to  find  an  expression  for  T*. 
Due  to  the  presence  of  the  factor  gj  in  all  these  expressions  and  also  the  equality  (2.66), 
we  can  find  the  constitutive  relation  for  However  before  doing  so  we  exploit  the  sym¬ 
metry  of  c'j*^  to  further  simplify  (3.3).  Since  c‘j“^  =  c‘j^  we  can  write 


cij^Kpa  =  ciif“K|to  =  c<i^  V  (3.20) 

Therefore, 


c‘j^Kp„  =  7  c‘j°«(Kap  +  ’fpo) 


(3.21) 


In  view  of  (3.21),  now  we  write  (3.3)  as 


T*  =  Tic/  J  J  g**'^  gj*  ^  d^  (3.22) 

With  the  explanation  presented  above,  now  we  write  the  constitutive  relation  for 


x«  =  {(1  +  ■^)n)c{j'^  +  .k  (1  +  ^)mc,^ 


3g 


^  Xi  (1  +  m  +  M  ±tS±2li)crrk/,^ 

4m  g  2m  '  3g  m 


+  (d  +  •^)  cjJ'“  +  y  3^  ('  +  4 
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+  -j-  —  (-ij  -  in)lc{l*‘)Kto 
8  g  m^ 


(323) 


If  we  further  restrict  ourselves  to  small  deformations  of  a  composite  with  initially  flat 
plies,  equation  (3.23)  can  be  simplified  to 


fc  2 

xy  =  {rncii^  +  (l-m)cii>^)7t/  +  ^  {mcP^  +  c^)K,a  (324) 

z  m 

Of  course  in  the  above  equations  no  distinction  should  be  made  between  covariant  and  contra- 
variant  components  of  tensors.  Using  (2.114),  equation  (3.24)  can  be  written  in  terms  of  dis¬ 
placement  vector  u  and  director  displacement  vector  S,  hence 


K  2 

Xij  =  {mCiSU  +  (l-m)^^}^!^/  +  y  {mcijii  +  qjglKto 

=  {mCiSop  +  (l-m)CiJ]j)  -j  (Uo^p  +  upa)  -K  2{mci^  +  (l-m)Ci^}(u3.o  +  6„)/2 

+  {mcg^  +  (l-m)Ci55^)53  +  y  (mCij^  +  Ci5S^}5/.o  (3-25) 

Using  the  symmetry  of  Cjjij/’s  (3.25)  can  be  written  as 

+  (l-m)Ci^p}Ua.p  +  {mcija3  +  (l-m)Ci55i^}u3,„ 

+  {mcijj^  +  (l-m)Ci^)6a  +  {naCj^^i  +  (l-m)CiS3U53 

+  Y  (mCiSU  +  c^^}5/.o 


=  {mCjjSj  +  (l-ni)CigJt)u,^o  +  {mciji^  +  (l-m)Ci5S)8k 


+  Y(>nCiii>  +  -ti^o®)5,« 
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The  same  steps  can  be  followed  to  calculate  S°^  and  we  record  the  procedure  here  KS 
S®  =  X  J  ^  T*“d4  =  --  I  g/d^ 

1  ^  1  ^ 

=  -f  J  48""  *i‘  d?  =  -f  J  4  +  2c-^fe 

^2  0  ^2  0 

+  c“^^Y33)gi*d^  =  Yk/  -p  J  +  K/p  J  g*''2c“i^g/d^  (3.27) 

This  is  basically  the  same  result  as  (3.22)  except  that  i  has  been  replaced  by  a  and  the 
integrand  has  been  multiplied  by  4-  The  first  integral  can  be  reduced  to  the  following  by 
refeiring  to  (3.18) 

gi/2  g  f  ^  [mcfj’^'  +  (~  -  m)  cf^]  +  [mcfj*^  +  (-^  -  m)cfj‘^] 

^2  m  og  tar 

+  m)^]  +  -^  >i [mcf'"  +  (-ij  -  rajcf"] )  (3.28) 

J  m^  og  m^ 

Similarly  the  second  integral  is  simplified  and  by  reference  to  (3.19)  the  result  is 


g^^  gj{%  [mcf^  +  (-^  -  m)c|‘j^]  +  [mcfj^  +  (-^  -  m)cfj^] 

3  nr^  og  m-’ 

1 3  1 4 

+  >4  [mcf^  +  (-\  -  m)cr^]  +  [tncf^  +  ( A  ' 

^  mr  i^g  m 

=  if  g''^  gjl(  J  +  +  >^i<7  + 


1  ,  1-m^ ,  .  1-m^  ,^ojip  .  ,  4 1  . 1-m^  .  l-m^ 


+  [*“  (— - )  + 

^3  ^  m2  ^  8g 


m*' 


4m2  lOg  m^ 


)C2®^)  (3.29) 
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Again  it  is  seen  that  because  of  the  factor  gj  and  the  relation  S“  =  S“j  gj  we  can 
readily  write  the  constitutive  relation  for  the  result  is 

s“'  =  +  -^)nicf>‘' 

2  m  3g  og 


This  is  the  general  constitutive  relation  for  in  linear  elasticity.  If,  as  before,  we  coniine 
ourselves  to  small  deformations  of  a  composite  with  initially  flat  plies  (3.30)  can  be  simplified 
to 


S"*  =  i-  (mefi"  +  ^  C2“J“)7u  +  4  ifCmc?**  +  ^  c?«>)K,  (3.31) 

2  m3  m-^ 

with  no  distinction  between  covariant  and  contravariant  tensors.  Once  written  in  terms  of  dis¬ 
placement  vector  and  director  displacement  and  simplified  as  done  in  obtaining  (3.26)  we  get 

S«i  =  i  cr»^)utf  *  j  c?i“)5k 

J  m  z  m 

+  4  5f(incf«>  +  cl-^bStp  (3.32) 

^  m* 

As  the  results  of  this  section  indicate,  even  in  the  simplest  cases  of  small  deformations  of 
an  initially  flat  composite  (composites  with  flat  plies)  higher  gradients  of  displacement  vector 
become  significant  and  they  appear  in  the  constitutive  relations  for  composite  stress  and  com¬ 
posite  stress  couple.  As  defined  by  (3. IS)  m  is  of  the  order  of  unity,  and  ^2  ^  usually 
small  lengths;  however  their  products  with  components  of  Cjju  and  even  the  product  of  their 
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higher  powers  with  elastic  constants  may  be  indeed  significant  quantities,  in  which  case  we 
get  ccmtrilHitions  to  Xy  and  Saj-  In  the  trivial  case  m  =  1,  =  ^2  0  we  get 

'*’ij  “  ®ijk/yic/*  Sotj  -  0  and  the  equations  of  linear  elasticity  are  recovered. 


BASE 
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4.0  COMPLETE  THEORY  FOR  LINEAR  ELASTIC  COMPOSITE  LAMINATES 

The  results  of  sections  (2)  and  (3)  are  combined  to  obtain  the  complete  equations  for  a 
linear  elastic  composite  laminate.  However,  before  doing  so  we  should  derive  appropriate 
expression  for  p^,  and  z^.  As  before,  we  assume  that  the  representative  micro-smicture  is 
composed  of  two  homogeneous  layers  with  respective  densities  Pi  and  P2  in  the  reference 
configuration  (p|  and  P2  are  constants).  Recalling  equations  (2.45)  and  (2.46)  we  write 

pg'« = p„G‘"  =  ;  pX'^d^  (4.1) 

^2  0 

Now  by  (2.13)  we  have 


G*1/2  =  g1/2(j  +  (42)* 

2G 

where  A  is  understood  to  be  the  sum  of  two  determinants  similar  to  those  expressed  in  (2.14) 
except  for  substituting  gjj  by  Gy.  We  have 


.  r  Pi  0<^<^i 

'’->■1  P2 

Substituting  (4.2)  and  (4.3)  in  (4.1)  and  using  (3.7)  we  set 


(4.3) 
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pg‘'V  =  p.G>'V  =  -^ 


(4.6) 


After  substituting  from  (4.2)  and  (4.3)  in  (4.5)  and  (4.6)  and  using  (3.17)  we  get 


PoZ^ 


r,,  ,  1-m^  ,  1-m^ ,  , 


3|,A 


Po^  =  —  ((1  +  +  ( 


l-m’  .  HiA  l-m< 


8G 


m-' 


)P2] 


(4.7) 


(4.8) 


For  a  composite  with  initially  flat  plies  (4.4),  (4.7)  and  (4.8)  are  reduced  respectively  to 


po  =  mp,  +  (l-m)p2  (4.9) 

poZ^  =  (mp,  +  P2)  (4.10) 

2  m 

PoZ^  *  %  (tnpi  +  P2)  (4*1 1) 

j  m 

To  formulate  the  complete  theory  it  is  also  worthwhile  to  derive  a  relation  between  the 
director  displacement  5  and  the  gradient  of  displacement  vector  u  in  the  6^-diiection.  In 
order  to  derive  such  a  relation  we  enforce  the  continuity  of  the  position  vectors  p*  and  P* 
between  two  adjacent  micro-structures.  Recalling  (2.1)  and  (2.2),  we  have  the  following  rela¬ 
tions  for  the  k***  micro-structure: 

p*(0“,03<*'\^,t)  =  r(e“.0^\t)  +  ^d(0“,e3(^),t)  (4.12) 

p*(0a  e3(k)  ^)  ^  R(ea,e3»))  +  ^(e“,e3»))  (4.13) 

Now  in  order  that  position  vectors  p*  and  P*  be  continuous  on  the  common  surface  between 
k*  and  (k+1)**  micro-siructures  we  should  have 

P*(e“,0^*'>,^2)  =  P*(9“.e^''*\0)  (4.14) 
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P*(0“.03(«44)  =  p*(0“.0^<‘‘-^^\O,t) 
Using  (4.12)  and  (4.13)  we  can  write  (4.14)  and  (4.15)  as 


(4.15) 


R(0“.03<>‘-^i>)  =  R(0“,03O'>)  +  ^2D(e“,0^‘‘^) 


(4.16) 


r(0“.03<>'+‘),t)  =  r(0“.03<^>,t)  + 


(4.17) 


Recalling  (2.98)  and  (2.99)  and  identifying  £u  and  eS  with  u  and  5  as  before  we  conclude 
from  (4.16)  and  (4.17)  the  following 


u(0“,03<k+‘).t)  =  u(0“.0^<^>.t)  + 


(4.18) 


5(0“  03(^\t)  =  {u(0“,0^^+‘\t)  -  u(0“,0^\t)] 

S2 


(4.19) 


By  smoothing  assumptions  and  approximating  the  right-hand  side  of  equation  (4.19)  as  the 
gradient  of  the  displacement  vector  in  the  0^  direction  we  have 


J(ga  03 ,)  ^ 

30^ 


(4.20) 


In  component  form  we  have 


6  =  5»gj  =  (ujgj)3  =  uj,3gj 


(4.21) 


6>  =  uj,3  ,  5j  =  Uj|3 

For  a  composite  with  initially  flat  plies  the  equation  (4.22)  reduces  to 


(4.22) 


8j  =  Up 


(4.23) 
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With  this  siinplification,  equations  (2.119)  reduce  to 

Yij  =  -J  (Ui.j  +  Uj.i)  (4.24) 

Kjo  =  Uj,ot3  (4.25) 

Using  (4.22),  equations  (2.121)  and  (2.125)  are  also  reduced  to 

Po  =  =  P(l+»i^  I  j)  (4.26) 

p  =  po(l-uj|j)  (4.27) 

The  constitutive  relations  (3.26)  and  (3.32)  for  ty  and  S^j  for  a  flat  composite  are  also  further 
simplified  by  using  (4.23) 

ty  =  {mCijil  +  (l-m)Ci(8)u,^/  +  {mCj^^^  +  Cj^)  Uk.a3  (4-28) 

and 

Saj  =  J  Ci^}uu  +  y  {niCajkp  +  cgplu^pj  (4.29) 

Now  we  can  write  the  field  equations  (2.113)  and  (2.124)  in  terms  of  the  displacement  vector 
u  and  its  gradients.  It  should  be  recalled  that  the  resulting  equations  are  the  linearized  field 
equation  for  small  deformations  of  a  composite  with  initially  flat  plies.  These  are  the  counter¬ 
pan  of  the  Navier-Cauchy  equations  in  linear  elasticity.  The  appropriate  equations  for  a  gen¬ 
eral  composite  will  be  derived  in  later  chapters.  Using  (4.9)-(4.11),  (4.23),  (4.28)  and  (4.29) 
we  write  (2.123)  and  (2.124)  as 
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{m4j^  +  (l-m)c4|[/luic.ta  +  {racijip  +  4]^p}  y  Uk,p3„ 

+  aj,3  +  [mpi  +  (l-in)p2]bj  =  [mpi  +  (l-m)p2]Uj 

+  ^  P2)Ui3 

and 

\  ^llncajlli  +  c^iS/IUtta  +  -y  l"Ka&  +  cij|!p)lltp3„ 

+  -  (OKiJi  +  (l-m)C®)UM  -  -I  (HKlii.  +  ^ 

1  l“"in^ 

+  [mpi  +  (l-m)pJCj  =  —  ^,(mpi  +  — y  p2)Uj 

+  ~  ^f(mpi  +  p2)up  (4.31) 

At  this  point  we  may  notice  that  an  ordinary  continuum  (a  single  material  continuum)  can  be 
regarded  as  the  limiting  case  of  a  composite  laminate  when  =  ^2  Therefore,  we  may 
anticipate  to  derive  the  equations  of  linear  elasticity  by  letting  m  =  1  and  ->  0  in  equations 
(4.30)  and  (4.31).  Doing  so,  equation  (4.30)  reduces  to 

Cojk/Ukjo  +  aj3  +  pbj  =  puj  (4.32) 

where  subscript  and  superscript  1  are  dropped  because  we  have  only  one  material.  To  sim¬ 
plify  equation  (4.31),  first  we  recall  the  definition  of  c  in  equation  (2.44)  and  notice  that  by 
the  mean-value  theorem,  c  -»  0  as  ^2  0,  hence 

Oj  -  C3jyUfc,;  =  0  (4.33) 
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Sul^tituting  for  Oj  fircmi  (4.33)  in  (4.32)  we  get 

Cajk/Uk.ta  +  C3jkiUij3  +  pbj  =  piij  (4.34) 

and  combining  the  first  and  the  second  terms  we  get 

Cijk/Uk.r,  +  Pbj  =  piij  (4.35) 

For  a  completely  isotropic  continuum 

Cijk/  ~  +  M-(5yt8j/  +  8i/Sji()  (4.36) 

where  X  and  p-  toe  Lame  constants  and  (4.35)  reduces  to: 


PUijj  +  (^■Hi)Uj.ij  +  Pbi  =  piii 

which  are  the  equations  of  motion  for  an  isotropic  media. 


(4.37) 


For  the  case  of  the  composite  laminate  we  can  also  eliminate  Oj  between  equations  (4.30) 
and  (4.31)  to  obtain  the  appropriate  equations  for  displacement  vector  u.  First  we  do  this  for 
a  static  problem  with  no  body  force.  For  such  a  case  we  let 


b  =  c  =  u’=  0 


(4.38) 


in  equations  (4.30)  and  (4.31),  hence 


i^aii  +  (l-m)cgi[,)Uk,ta  +  (mc^j 


1-m^ 

m 


Ca?p)  'Y  =  0 


and 
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■y  (mcojili  +  Coji[/}Uk.ta  +  ■—  {m4ji[p  +  45[p}Uk,p3a 

^  in  j 

K  2 

+  aj  -  {mcij)/  +  (l-m)cg)/)uk./  -  -j  ciji?o)uk.a3  =  0 

Eliminating  Oj  between  these  equations,  we  get 

{mc<ji[i  +  (l-m)c®/)Uk,ta  +  {mc^jilp  +  4^p)  Y  “upa 

+  {mc^Jil,  +  (l-m)c^J^/}uk,o  +  —  ci51o}Uk,a33 

-  Y  {m4jl[/  +  cg^/)Uk./a3  -  Y  cSp)Uk,pa33  =  0  (4.39) 

By  combining  the  first  and  third  and  also  the  second  and  fourth  terms  of  the  equation  (4.39), 
we  get 

t  2 

{mcjjljl  +  (l~m)Ci(S)utfi  +  Y  {mcijl,  +  cg>a)Uk.ia3 

-  Y  C§/)Uk,fa3  -  Y  4ji[p)Uk.pa33  =  0  (4-40) 

This  is  a  fourth  order  partial  differential  equation  for  displacement  vector  u.  Now  we  apply 
this  equation  to  a  composite  laminate  whose  micro-structure  is  composed  of  two  isotropic 
layers.  For  such  a  case  we  can  write 

Cijkl  =  ^iSij5ic/  +  lAi(5ik5j/  +  5i/5jjt)  (4.41) 

Ci^)  =  +  M.2(5ik5j/  +  Sj/djk)  (4.42) 

where  X4  and  m  (i  =  1,2)  are  Lame’s  constants  for  the  respective  layers.  Introducing  equa¬ 
tions  (4.41)  and  (4.42)  in  (4.40)  we  get  the  following  for  each  term: 


{inc®  +  =  m{XiUtkj  +  ^i(utjk  +  Uj.a)} 

+  (l-m){X2Uijy  +  ^2(“k.jk  +  Uj.//)} 


=  {in(Xj  +  M-i)  +  (l-m)(X2+42)}“k.kj  +  +  (l“n^)^2)Uj,// 

^  ^ijkoJ^kjotS  ~  +  M^i(S^5jQ  +  Sia^jk))^k.i(x3 

+  (MijSto  +  42(5ik8ja  +  8ia5jk))Uua3 


-  ni{X,i5ij6potUp  iot3}  +  mni{5iij8jaUij^ioi3  +  Sp^SjijU^  po3 } 


+  iLSlI  X.2(5ij5paUp,ia3)  +  M5|k8jaUk.io3  +  SpoSjkUi^ios) 


m 


{mXi  +  ^2}Vjo3  +  {»"^l  +  l^2}8jaUk,kot3 


m 


m 


+  {m^i  +  '  H2}Uj,aa3 


^  ~  n^{X,i5otj5i{/  +  M^i(5Qjt5j/  +  5a/8jk)uk,/a3 


+  ^  {^28aj8k/  +  ^2(8oac8j/  +  8„^jlj))Uk,/a3 


-  mA,i(8ajUk,ko3)  +  nTM.i(Uaja3  +  Ujaos) 

.  (1-m^)  ^  X  .  (1-m^) 

;;;  ^(8ojUk,ka3)  +  — -  MUajaS  +  Uj  003) 

=  (mXi  +  ^  ^  ^2)8ajUk.ka3  +  ^l2)(Ua^ja3  +  Uj  „a3) 

m  ■’  m  ■’ 


(4.43) 


(4.44) 


(4.45) 


BAS 


(mCojip 


1-m^ 


^o^)“k,pa33  ~  ■*■  +  ^opSjk))^lc.Pa33 


1-m^ 


m“ 


2~  {^Saj^kp  +  M^lC^akSjp  +  5ap8jk))Uk.poi33 


=  mXiSajUp  ptt33  +  n4il(SjpUa,pa33  +  Uj.ao33) 

+  ^2  ^^aj“p.po33  +  ^  M-2(Sjp“ot.Pot33  +  “i.aa33) 

=  (mX,i  +  T-  ^)SttjUp  p<jt33  +  (mUi  +  —  |l2)(5jpUa,pa33  +  Uj^3)  (4.46) 

m  m 

Substituting  (4.43)-(4.46)  in  (4.40)  we  obtain 

-j-  {mXi  +  ^  ^  ^  ^)Uaja3  +  {™4l  +  "  i^}(5joUkJta3  +  “j.ao3) 

+  (m(Xi  +  Ji,)  +  (l-m)(X2  +  H2))uiaj  +  +  (l-m)^2}uy/ 

-  Y  (mXi  +  — ^8ajUkj,a3  -  Y  (mui  +  H2)(Ua.ja3  +  Uj.ao3) 

-  Y  ("^1  ™  ^)8ojUp.pa33 

m 

%\  1  ^ 

-  Y  ■*■  — T~  ^^2)(8jpUo,pot33  +  Uj  aa33)  =  0 

j  m 


Now  we  introduce  the  following  definitions  in  the  last  equation; 
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X,i2  =  mXi  +  (l-ni)A2 


Hi2  =  m|i,  +  (l-m)H2 


Xi2  -  niA.1  +  ^  ™  X2 
m 


.  l-*n^ 

1X12  =  mUi  +  — — 
tn 


(4.47) 


Xj2  —  niX-i  + - —  X2 


m* 


)ii2  =  mui  + 


1-m" 


M2 


The  result  will  be 


— 

(Xi2  +  lti2)Uic,icj  +  +  "Y  {^12“oga3  ^  Pl2(SjaUkJta3  +  “j.ao3)) 

“  "Y  ^12  SojUk.ka3  “  “Y  ^12(“a,jo3  +  Uj.ao3) 

- 

- T  ^12SojUp.Po33 - T  i^l2(SjpUo,po33  +  Uj,ao33  =  ^  (4.48) 


The  above  equations  are  counterpans  of  the  classical  equations  for  linear  elasto-static  prob¬ 
lems  in  the  absence  of  body  forces.  We  need  to  write  these  equations  in  the  expanded  form. 
The  result  would  be  three  scalar  equations  as  follows: 
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^  a  dui  3u2  9u2  9^1  9^ut  9\ii 

(Xi2  +  1^12)  ,  (  -H  *  -f- - r)  +  (li2(  '  "  /,  +  '  "T  +  ■  «  ) 

P12/ V9ei  ^2  303^  ^‘2^302  002  dej^ 


.  Il  r  9^  ^Ui  ^“2  _  9^  ^Ui  9u2  9u3 

2  90,903  '^  902^'^  2  90,903  902  903^ 


9  9^ui  9^u,  4i  z-  9^  3ui  9u2  9u3 

T  +  -^)  -  Y  ^12  30j^  (-^  ■*■  907 


^1  _  9^  .  3ui  ^  ^“2 ,  ^1  -  _9_  /  ^1  ^  , 

2  ae^ae^  ^001  002^  2  003  ^002 


_  1l  f  /in  in ^  _  1l  =  /in  in\ 

3  00,003  90,  902  3  90,903 


i2u. 


=  jj_  ,iiii_  ^  . 

3  0032  ^  002  0022^ 


(4.49) 


This  is  the  first  equation  for  j  =  1  which  can  further  be  simplified  ar> 


^  dUl  '^*•'7  vwa 

(X,i2  +  n,2)  +  -^)  -»-  ^i2(ttt  +  “rrr  +  ttt) 


9u2 


9u3 

90^ 


9^,  9^,  9^, 

Bdi  902  90-? 


4,  —  9^U3 


/T  ^  ^  9^  ,  3u,  9u2 

3  ^  36,3032  <ae,  363' 


4?  =  92  „ 

3  3632  *  36,2  36,2  * 


(4.50) 


The  second  equation  will  be 


(Xi2  +  »ii2)  +  ■^)  +  ^1 

0X92  ^^3 


^1  _  —  ^^“3 

T  aa^ 

^  OOoOXj'i 


^  =  .  3^  .^“1 


=  a^  ^2  ^ 

3  afl2  afl2  ^  aa^  ^ 

J  OO3  OD|  OD2 


And  the  third  equation  is 


,  a  3ui  9u2  du3 

(X,J  +  n,2)  -jj-C-jg-  +  +  a^-)  +  •‘n 
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a  dui  3u2  dua  d^3  9^U3  d^3 

(Xi2  +  ^12)  (■5S“  ■*■  ■‘'  ■5S")  ^  + 


363  301  302  303 


30j^  302  303 


^1  _  32  3ui  3u2 


=  32  3^U3  ^^3 .  _  _ 

3  ae^  38?  3? 


(4.52) 


For  future  reference  we  will  also  calculate  Oj  for  static  problem  in  the  absence  of  body 
forces.  The  preceding  results  are  used  in  conjunction  with  equation  (4.31)  which  has  been 
simplified  for  such  a  case  and  reproduced  before  the  equation  (4.39).  From  (4.45)  and  (4.46) 
and  (4.47)  we  have 


{mciji|/  + 


(l~m^) 

m 


Co]ici)Uk,ta  =  ^12  SojUjuka  +  Pl2(Uajo  +  ‘^j.ota) 


(4.53) 


^ojip)’J'k,33o  “  ^12  ^oj‘*p,pa3  M'12(^jp‘*o,pa3  '^j.oixs) 


(4.54) 


Using  (4.36)  we  get 


{mcj^  +  (l-m)cjjii}Ukj  =  [mXj  +  (l-m)X2]53j5,j^;^  / 


+  [m^i  +  (l-m)p2][83k5j/  +  53i5jk]Uk,i  =  Xi253jUy^  +  Pi2(u3j  +  up)  (4.55) 


cijico)Uk,a3  -  ^1253j5ta“k.a3  +  M^12(53k5ja  +  53aSjk)Uk.a3 


=  Xi263jUa,a3  +  12l25jaU3.a3  (4.56) 

Substituting  (4.53)-(4.56)  in  relation  for  Oj  we  obtain 
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*  ^lAjUkOc  +  ^ll2(“3o  +  UjJ)  +  ”^[^1283jUa.o3  +  Pu^joUS^is] 


^  f^^aj^kjta  ■*■  J^12(*Vija  "*■  “j,aa)l 


3 


- [^12SajUp.pa3  +  ^12(SjpUa,po3  +  Uj^oos)! 


(4.57) 


Writing  down  the  components  of  Oj  separately  we  get 


_  ^  du3  9ui  ^  4i  „  ^3  5r  9  ^“2  ^“3  ^ 

^  ae7^  T  ael^  "  T  ^ 


^1  ^  a  dui  ^  9u2,  ^1  ^  ,d^ui  a% 

2  301  ^301  302  ^  2  ^^^^30f  302^ 


^?  =  32  ,  ^1  .  3u2 


3  301303  ^301  302 


+  '^) 


11  =  ^  .^“1  .  in.  _  11  =  -i- 

3  301303  ^301  302^  3  303  ^902  002^ 

3ui  3u3  &i  _  3^U3  ^1  —  3  5ui  3u2 

1^12  ("^  +  ■^)  +  “T  (1^12  ~  ^12)  T  (^12  +  M^12)  +  ■^) 

CrU3  dvj  Z  CWjCWj  Z  C^l  (/V| 


^1  _  ,  ^^Ui  4?  -  =  ,  32 


3ui  3u2 

30i303 


4?  3  ,32Ui  32U1-= 

3  363  30?  ■"  301 


(4.58) 


3u2  ,  3u3  ,  4i  .  3  ,  ^“1  .  ^“2 


O2  =  ^^12('^  +  ■^)  +  ^  (^^12  “  ^12)  no  na 
Ov3  Ov2  Z  C/\/2^^3 


-  ^  (^12  + 1^12)  (■^  +  ■^) 

Z  Ct72  *2  CFO2 
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^  ^2.  ^  au2^ 

■  T  ^  -  T  ^  a^  aeT^ 


_1l  ^  3^2  ^  ^U2  - 

3  ae,  ^  302  302  ^^"12 


(4.59) 


^1  4i  4i  — 

O3  =  Xi2Uk,k  +  2M^i2U3j  +  —  (^12“o.o3)  “  Pl2(“ojo  +  “3,00) - ^  )2l2“3.ao3 


_  aui  au2  au3  au3  3  aui  au2 

“  35“  ^^^12  '55“  ■•■  ^12  (ts”  +  ”55“) 

CrU|  (7U2  CKI3  C^3  Z  (/D3  aU|  (7U2 


4i  3  aul  au2  a^3  a^u3 


-I- 


a\i3 


,  aui  au2  ^  au3  ^1  :t-  a  ^1  ^“2 

®3  ®  ^i2('55"  ■’■  "55“)  ■•■  (^12  +  2412)  -55“  ■•'  “5"  (^12  -  F12)  *55“  (■55"  ■*•  ■55“ 

<^1  CrU^  Z  <7U^  vUi  0XJ7 
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4i  a\i3  a^u3  _  3  a^3  a^u3 

-  T V  "  V  ■  T  ^  <liF  *  liF> 


(4.60) 


The  constitutive  relations  (4.28)  and  (4.29)  for  x^j  and  S^j  for  a  flat  composite  are  also 
simplified  here  for  a  composite  laminate  whose  micro-structure  is  composed  of  two  isotropic 
layers.  Using  (4.36)  first  we  simplify  different  terms  of  the  expressions  (4.28)  and  (4.29) 


{mcijj^-i*  (l-m)cg)}U|t,,=  {n^i  +  (l-m)A^}5ij5,^/u,^,/ 

+  {mu,  +  (l-m)^2)(5ik5J/  -i-  fij/fij,  +  u^,/ 

=  ^i25ijUk.k  +  H,2(Uij  +  uj,i)  (4.61) 


BASE 


wtere  (30.47)|^  are  used. 


2 

{mcijil;  +  — =  ^nSajSkiUu  +  Wni^ak^jl  +  ^afiik)^U 

=  ^125ajUkJc  +  Pl2(Ua,j  +  Uj.a) 

2 

{mCiJit  +  ^  cg>,)Ui^  =  ^uSijSka'ik.etS  +  Pl2(SikSja  +  MjkKoS 

=  ^i25ijUa,o3  ■*■  Pl2(^ja*^i,a3  ^iaUj,a3) 

l_Uj3  _  _ 

(mciiJ,  +  =  i, Aj8kpu^3,  +  ii,,(S<^8ip  + 

=  Xi28ojUp.p3  +  Pl2(8jpUa.p3  +  Uj.o3) 
Substituting  (4.61)-(4.64)  in  (4.28)  and  (4.29)  we  obtain 


Xij  =  XijdijUy,  +  ^i2(Uij  +  Uj,i) 

+  ~2  I^12Sij“a.o3  +  Pl2(SjaUi,a3  +  ^ioUj^oa)! 

— 

Soj  =  “  [^12^aj'*k.k  ■*■  Pi  2(^00  ■•■ 


(4.62) 


(4.63) 


(4.64) 


(4.65) 


+  "J  %\  l^l2Saj“p,p3  +  Pl2(SjP“o,P3  +  '^j,a3)J 


(4.66) 
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5.0  LINEAR  CONSTITUTIVE  RELATIONS  FOR  A  MULTI-CONSTITUENT 
COMPOSITE 

In  this  section  we  assume  that  the  representative  micro-structure  is  composed  of  n  layers 
with  different  constituents.  For  such  a  micro-structure  we  let 

VcS  =  (a  =  U . n)  (5.1) 

where  C(^  (a  =  l,...,n)  are  material  constants  in  the  associated  layers.  As  before  the  variable 
^  is  designated  to  change  across  the  micro-structure  whose  thickness  is  assumed  to  be  It 
should  be  noted  that  although  the  micro-structure  is  composed  of  n  layers,  ^  is  still  supposed 
to  be  a  very  small  number.  The  range  of  variation  of  ^  in  the  /***  layer  of  the  microstructure  is 

from  ^^1  to  4/  where  /  =  l,...,n  and  ^  =  0.  This  convention  is  adopted  due  to  its  agreement 

with  the  special  case  of  a  two-layered  micro-structure  which  was  studied  before.  We  further 
define  (n-1)  constants  mi, . . .  ,nvi  according  to  the  following  relations 

=  m^,  (/  =  l,...,n-l)  (5.2) 

As  a  result  of  this  definition  the  thickness  of  the  i*  layer  of  the  micro-structure  is  equal  to 
(m/-m/_i)^  where  /  =  l,...,n  and  m^  =  0,  m„  =  1.  The  composite  stress  vector  T'  and  the 
composite  couple  stress  S®  and  other  quantities  are  obviously  defined  over  the  whole  thick¬ 
ness  of  the  micro-structure.  For  example  we  have 

I  ^ 

T*  =  J  T**d^  (5.3) 

Sn  0 
1  ^ 

S“  =  J  ^  T*“d^  (5.4) 

Sn  0 

In  order  to  derive  appropriate  constitutive  relations  we  make  the  following  definition. 
Let 


BASE 
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a  =  a/  for  4i-i  <  ^  <  4/  (/  =  1.2,....n)  (5J) 

where  a/s  are  constants  and  ^  =  0  as  noted  earlier.  The  function  a  is  piecewise  continuous 

for  4  €  (0,  and  we  can  evaluate  the  following  integral 

Using  definitions  (5.2)  we  further  simplify  (5.6) 

J  ^^ad^  =  -j~  2  a^m/‘'-^'-m]iVl^n^‘"*  =  -j—  S  a/(m/''+’-m/‘iV)  (k  -1)  (5.7) 

where  nio  =  0  and  nin  =  *^0  simplify  the  final  results  in  constitutive  relations  we  first  notice 

that  the  integrals  which  appear  in  these  equations  are  the  weighted  averages  of  the  constitutive 
coefficients.  So  we  adopt  the  following  definition 


1  ^ 

l(k)pqrsA  J  ^kgpqiSjj^ 

^  0 


which  by  (5.7)  is  seen  to  be  equal  to 


(5.8) 


2  -  m,'!:!')  (5.9) 

K »  1 

We  use  the  same  contravariant  or  covariant  index  notations  for  I  and  c.  However,  the  weight¬ 
ing  number  k  is  always  written  as  a  superscript  in  parentheses.  Whenever  the  covariant  com¬ 
ponents  of  constitutive  coefficients  are  used,  the  layer  index  (/)  is  also  wrinen  as  a  superscript 
in  parentheses.  Recall  (3.22)  which  for  the  present  situation  can  be  written  as 

T*  =  Yk/-p  J  gj*  d^  +  Kta-p  J  (5.10) 

Combining  (3.4)  and  (2.9)3  obtain 
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gp  =  (1  +  Ij^gp  +  ^^gj) 


(5.11) 


and  (3.5)  reads 


g3  =  g^^  (1  +  •|j)g3 

The  first  integral  in  relation  (5.10)  is  simplified  using  (5.11)  and  (5.12).  We  have 


(5.12) 


1  J  g-1/2  g.-  ±  J  (ciUki  g-1/2  Jg-  H.  ,i3U  g.1/2 

WO  ^0 


^  ^  +  g^'^  >-pgj  J^d  +  ■^)c‘^‘^d^ 


+  g*'^  g3  J  (1  + 


2g 


"  (5.13) 

Now  we  use  definition  (5.8)  to  simplify  (5.13).  The  following  expression  would  be  the  result 


1“  term  in  (5.10)  =  +  —  l»)iPW)j  (5  ^4) 

2g  2g 

The  second  integral  in  (5.10)  is  also  simplified  similarly 


1  ^  ■  1  ^ 

—  j  ^c'^“  g**'^  g/  =  T-  J  ^  g*''^  (c'^**gp*  +  c^**g3)d^ 

^  0  ^0 


=  -^  I  4  g'**  (1  +  ■"  I  27’*^ 


+  /  4  g''^  (1  +  ■^)c“'“g3<i4 
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!U  4. 

=  g''^  gjt-f  I  50  +  47X=‘^‘15  +  ^  -f  J  5^0  +  •l7)c‘^“<i5) 

Sn  0  Sn  0 


(5.15) 


which  again  by  using  (5.8)  is  simplified  to 


2™*  term  in  (5.10)  =  ^  (5.16) 

J  2g  P  2g 

Substituting  (5.14)  and  (5.16)  in  (5.10)  we  get  the  following  constitutive  relation  for  T* 


T  =  +  A  i(i)ijw  +  ;^j(i(i)iPk/  +  A  i(2)iPk/)jY^ 

+  [l(l)ijto  +  A  l{2)ijto  +  gj  (5.17) 

2g  2g 

The  expression  inside  the  bracket  is  obviously  the  constitutive  relation  for  I'j.  Hence 


xij  =  ii<o)ijw  4.  A  iwm  + 


+  [l(l)ijfa  +  A  l(2)ij/a  +  xj(l(2)ipto  +  A  i(3)ipta)j^^^ 


(5.18) 


The  same  steps  are  followed  to  derive  the  constitutive  relation  for  the  composite  couple 
stress  S“.  By  (3.27)  we  have 


S“  =  Yk,  A I  ^  c“J‘'^g*d^  '^/p  "r  J 

Sn  0  Sn  0 

which  can  be  reduced  to  the  following  form  by  exactly  using  the  same  procedure 


S“  =  {lO)“jW  +  A  l(2)ajk/  ^  jj^j(2)a7k/  ^  A  l(3)onfk/)]Y 

2g  ^  2g 


+  (I(2)aj/p  ^  A  i(3)aj/p  ^  ^iO)ayl^  ^  A  l(4)a#)]K^}  gl/2  g. 


(5.20) 


BASE 


Subsequently  the  constitutive  relation  for  would  be 


S“j  =  +  A-  i(2)ajk/  +  j^^j(2)a‘)(k/  ^  j(3)a'>)c/jjY^ 

2g  2g 


+  [i(2)aj/p  ^  A  l(3)oj/p  ;^(i(3)ccy^  +  A  i(4Xryip)j^^p 


(5.21) 


For  small  deformations  of  a  composite  with  initially  flat  plies,  the  foregoing  equations 
are  further  simplified.  The  resulting  relations  are  recorded  here: 


(5.22) 

S“j  =  (5.23) 

with  no  distinction  between  contravariant  and  covariant  components.  In  terms  of  displace¬ 
ment  vector  u  and  its  gradients  these  equations  can  be  written  as  follows: 


=  ^ijkW/  'ijlaU/.a3  (5.24) 

^aj  “  ^aj^^k./  ^aj^P^/.P3  (5.25) 

Using  (5.9)  the  constitutive  coefficients  are  written  in  the  expanded  form 


Ci^nV  -  nV-i) 

I®  =  7  c5l>y(m2  -  nvl,)  (5.26) 

=  7  ciii(nv^  -  nvli) 

Z  1=1 


To  recapitulate,  d^^ih  (r  =  l,...,n)  are  the  constitutive  coefficients  of  the  micro-structure  layers 

and  nij’s  (r  =  1 . n-1)  are  dimensionless  constants  related  to  the  thicknesses  of  different 

layers  with  nio  =  0  and  nin  =  1 . 
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If  the  micro-structure  is  composed  of  n  isotropic  layers,  we  can  write  (5.26)  in  terms  of 
the  Lame’s  constants  of  various  layers.  In  fact,  we  have 

Cijli  =  ^rAjSk/  +  M<r)(Sik5j/  +  Sygjij)  r  =  l,...,n  (5.27) 

For  such  a  case,  the  relations  (5.24)  and  (5.25)  are  written  in  expanded  forms  as  follows 


=  “kjcSy  2  >.(,)(ni,-nv_i)  +  (Ui^+Uj,;)  S  4(,)(mi-nv_,) 

r=l  fsrl 


E  n 

+  Y  ^^iAaU/,a3 


^  Hcd  2  Ii(r)(m2-m2.,)] 
r=l 


'^kje  Sy  J*  \r)  "*■  ” 

^  n 

+  ^  (Ui.a3Sja+Uj.a35ia)  ^  li(r)W  (5.28) 

1  "  -  n 

Saj  =  "T  4n[5oj5k/Uk,/  2  X(j)Anv  +  (5aic5j/  +  5Q^j|()uk  /  I  P(j)Amj^] 

r=l  1=1 

1  2  ”  -1  " 

+  T  [5ajS/pU/.p3  Z  X(^)AnV  +  (5ay5jp  +  5ap6j/)U/,p3  Z  li(r)Anrv^] 

^  rsl 


j  n  n 

Saj  =  7  ^[5ajU|c.k  2:  +  (Ua.j+Oj.a)  2:  ll(r)Ani2] 

^  r*!  isj 

i  n  n 

+  T  W[SajUp.p3  2:  A.(,)Am3  +  (5jpUo,p3  +  Uj,a3)  Z  li(,)Anv’] 

J  r=\  r=l 


(5.29) 
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where  for  brevity  we  have  introduced 

AmP  =  mP  -  tnP.1  r  =  (5.30) 

and  in  the  above  relations  p  =  1,2,3. 

In  order  to  obtain  the  complete  field  equations  for  a  linear  elastic  composite  whose 
micro-structures  comprise  n  layers  we  should  substitute  the  constitutive  relations  (5.24)  and 
(5.25)  in  the  equations  of  motion  (2.123)  and  (2.124).  However  before  that  we  should  obtain 
appropriate  expressions  for  p,,,  and  PqZ^.  We  assume  that  the  micro-structure  layers  are 
homogeneous  with  densities  p^^^  (r  =  l,...,n)  in  the  reference  configuration.  Recalling  (2.45) 
and  (2.46)  we  can  write 

pg^^  =  poG‘^  =  j-  j  (5.31) 

^  0 

where 

Po  “  Po^^  (or  ^1-1  <  ^  ^  (r  =  l....,n) 

and  (5.32) 

^  =  0 

Using  (4.2),  (5.32)  and  (5.7)  in  (5.31)  we  conclude 

p»  =  r  J  2  PfM  +  ^  P» (5.33) 

^  0  2G  1=1  4G  r=l 

We  will  proceed  similarly  to  calculate  p^z^  and  p^z^  using  their  definitions  for  the  present 
situation.  By  (2.42)  we  have 

,  ^ 

pgi'^z'  =  p„G''^z'  =  j-  j  ^poV'^d^  (5.34) 
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^  0 


Similar  to  what  was  done  in  the  derivation  of  (5.33)  we  write 


(5.35) 


TO  fe  2 

Po*'  “  ■  "2  ^  ^1  P“ ■*■  Ji  (5-3«) 


°  ^  4^(1  +  ■f|w<'4  =  1 4,^  I,  Po^'^W  +  -^  I,  pi'’M‘  (5-37) 


For  a  composite  with  initially  flat  plies,  equations  (5.35),  (5.36)  and  (5.37)  are  reduced 
respectively  to 


Po  =  S 
r=l 


PoZ^  =  T  Po'W 

I  r=l 


(5.38) 


(5.39) 


PoZ^  =  -7  \  Po' W 

j  1=1 


(5.40) 


Now  we  can  substitute  (5.24),  (5.25),  (5.38)-(5.40)  and  (4.23)  in  (2.123)  and  (2.124)  to 
derive  the  linear  equations  of  motion  for  a  flat  composite.  The  resulting  equations  are 
recorded  below: 


n 

Ioji[/Uk,ta  +  Ujiu/.ap3  +  bj  +  aj,3 

=  iij  Z  pW/tal,  +  i  4„up  Z  p<'>Am,7 

T=1  Z  ■’  I=l 


(5.41) 
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Io]^“k.ta  +  Ujiu/,ap3  +  Oj  -  p3 

+  Cj  I.  p \  S  +  1  ^2uj.3  Z  pi^^Am^  (5.42) 

The  appropriate  differential  equation  for  the  displacement  vector  u  is  obtained  by  elim¬ 
inating  Oj  from  equations  (5.41)  and  (5.42).  For  static  problems  with  no  body  force  these 
equations  reduce  to 


^aj^P^lc,ot3P  "*■  *^j.3  ^ 

^aji^P^k.ap3  +  Oj  ~  lijlci^ic,/  “  lijlcP*^kp3  “  0 
Eliminating  Oj  between  these  equations,  we  get 

+  Iaj^Uk,op3  +  lijk^/Uic./3  +  li^pu,j^p33  ”  Uj^/“k./a3  “  Iaj^pUk.aP33  =  ^ 
or 

'tijS“k.fi  +  *ijktUk.ip3  “  Iaj^/“k./a3  “  U]^pllk.aP33  =  0  (5.45) 

This  is  a  fourth  order  partial  differential  equation  for  the  displacement  vector  u.  The  constitu¬ 
tive  coefficients  (r  =  0,1,2)  are  already  written  in  expanded  forms  in  Eqs.  (5.26). 

For  a  composite  laminate  whose  micro-structure  is  composed  of  n  isotropic  layers  we  use 
(5.27)  and  (5.26)  to  rewrite  (5.45).  The  result  is 

n  n 

5ij8klUk./i  +  (SikSj/  +  5i/5jk)Uk.fi  P(r 

+  Y  ^(5ij5kpUk.ip3  ^(rjAnv^  +  (SjkSjp  +  5ip5jk)Uk,ip3  P(r)Anv  ] 

j  n  n 

-  T  ^l5„j6k,Uk  to3  ^  hA^  (^otkSj/  +  5oj5jk)Uk.ia3  2:  P(,)Am2] 

Z  T=\  r=l 


(5.43) 

(5.44) 
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_  1  J:2| 


1  ^  -I  ”  -I 

T  ^fSojSkpUtoP33  ^  ^(ry^  +  (Sok^jp  +  Sop5jij)Ui5^3  I  =  0 

1*1  Ifel 


Ui,ji(  Z^  X(r)Anv  +  Z^  m^jAirv)  +  Uj.ii  Z^ 


2  11  n 

3  4n[Upop3  ^r)Ani2  +  (5jpUyp3  +  Uj,pp3)  Z^  ^(r)^!!!^] 


2  II  II 

-  7  4n[SajUu:a3  ^  +  (Ua.ja3  +  Uj.ao3)  ^  ^(r)Am2] 

^  1*1  1*1 


-  T  ^n[5ajUp.ap33  2:  \r)^  +  (SjpUa,ap33  +  Uj.ota33) 

J  r*l  r=l 


“Mi  (^(r)  +  mr))Amr  +  Uj^jj  Z  H(r)Amr 
1=1  f=l 


1  " 

■*■  7  ^^^“aoa3  ~  ^aj'ii.ias)  ^  t^r)  “ 

^  I*l 

1  7  "  ” 

-  7  ^(SojUp,ap33  Z^  [\,y  +  ^(r)]Am?  +  Uj,ao33  =  0  (5-46) 

There  are  three  partial  differential  equations  of  fourth  order  for  displacement  vector  u  and  we 
can  write  them  in  two  separate  sets,  one  for  j  =  Y  the  other  for  j  =  3.  For  j  =  y  we  obtain 


n  n 

“i-)i  (^(T)  +  + ^.w  21 

1=1  1=1 

1  " 

“  7  4n  “3,733  (^r)  “  M^(r))^“™r^ 

It”  " 

“  7  ^^“P.7P33  (^(r)  +  +  “7.aa33  =  0  (5.47) 

3  r*l  r*l 

and  for  j  =  3  we  get 


BASE 


(5.48) 


“i.i3  (^(r)  +  )i(r))^  +  ^S.ii  ^  ^“o.a33  (\t)  “  ^i(r))^^ 


-  y  ^nU3,a„33  H(,)Am3  =  0 

For  the  special  case  where  n  =  2  these  equations  reduce  to  the  equations  (4.50)-(4.52)  derived 
before  for  a  bi-latninatc  composite.  From  (5.44)  we  can  also  calculate  Oj 


which  for  the  isotropic  case  reduces  to 


(5.49) 


••  II 

=  (SsjSjciUjt /)  1  ^(r)An\  +  (53|t5j/  +  Z  M-(r)^™r 

r=l  nsl 

“  -r  tSajSk/Uk./a  ^  +  (SakSj/  +  ^  li(rA"r^] 

*  f=l  1=1 


■*■  ”  t^3jSkpUk.p3  +  (^SpSjk  +  83kSjp)Uk.P3  1^(1 

“  *T  ^  lSoj5kpkk,aP3  ^  +  (5(jk5jp  +  5op5jk)Uk  „p3  Z  p^yjAnij^] 

J  r=l  1=1 

The  equations  (5.50)  are  written  in  two  separate  sets.  For  j  =  y  we  get 


(5.50) 


n  ^  n 

CJy  =  (lly,3  +  U3  y)  ll(r)^mj  +  —  U3  y3  -  X(j)]Amr^ 


2  t“o.a7  (^(r)  ■*■  +  Uy.aa  M^(r)^™r^] 


_1_ 

3 


1 

^  ['^0.073  (^(r)  1*7,003 


r=l 


(5.51) 


and  tor  j  =  3  we  obtain 


(5.52) 


03  =  Ujtjc  Z  +  2u3,3 

r®l  fsl 

-  %  ^“030  +  «3.ao)  2:  ^l(r>Am2  +  ^  Uot,o3  \t)^ 

^  f=l  Z  1=1 

1  " 

-  J  ^^3.aa3 

For  the  special  case  of  a  laminate  with  two  layer  micro-structure  (n=2)  these  equations  reduce 
to  (4.58)-(4.60). 
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6.0  DEVELOPMENT  OF  THERMO-MECHANICAL  THEORY  FOR  COMPOSITE 
LAMINATES 

In  order  to  develop  a  thermo-mechanical  theory  for  the  composite  laminates,  we  begin  by 
writing  down  the  local  balance  of  energy  and  the  Clausius-Duhem  inequality  for  the 
representative  micro-structure.  First  we  introduce  the  following  additional  five  quantities 
which  we  associate  with  a  motion  of  the  micro-structure: 

The  specific  internal  energy  e*  =  e*(0“,6^\^,t) 

The  heat  flux  vector  q*  =  q*(0“,0^^\4,t) 

The  heat  supply  or  heat  absorption  r*  =  r*(0“,0^^^^,t) 

The  specific  entropy  “n*  =  Ti*(0“,0^^‘‘\4,t)  and. 

The  local  temperature  0*  =  0*(0“,0^^\^,t)  which  is  assumed  to  be  always  positive.  The 
equation  for  the  local  balance  of  energy  —  the  first  law  of  thermodynamics  —  can  be  written 
in  the  following  form 


•  •  *ii'  *  ‘k  /\ 

p  r  -  p  e  +  X  %  -  q  Ik  =  0 

where  p*  is  the  density  of  the  micro-structure,  q**'  and  y-*  are  defined  by 


(6.1) 


q 


•k  * 

=  q  gk 


•  •  1  .  • 
Yij  “  2 


(6.2) 


and  covariant  differentiation  is  performed  with  repsect  to  the  metric  tensor  gj*  of  the  micro¬ 
structure.  Recalling  the  relations 


gi  =Vi 
gij  “  gi  '  gj 


(6.3) 


A  S  E 
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T*i  =  g-i/2c-ijg- 

We  can  write 


*•  !•*  1,*  *.  •  *. 
Yij  =  7  gij  =  7  (Vi  •  gj  +  Vj  •  gi ) 


2  2 

Using  (6.3)3  the  symmetry  of  x*'^,  we  can  write 


(6.4) 


=  T  (v.i  •  x*yg*  +  V*  •  x*yg*) 


=  T  (''.i  ■  +  V*  •  g*-i/2T*j)  =  g*-i/2T*i .  V ; 


(6.5) 


As  for  the  divergence  of  the  heat  flux  vector,  we  have 


div  q*  =  (g*t/2  q*k)^ 

g 

Introducing  the  results  (6.5)  and  (6.6)  in  (6.1),  we  can  write  the  local  energy  equation  in  the 
following  alternative  form 


p*r*  -  p*e*  +  g*  .  y  *  _  (g*i/2q*k)  _  q 

The  energy  equation  can  also  be  written  in  terms  of  the  Helmholtz  free  energy  function 
defined  by 


V*  =  e*  -  0*Ti*  (6.8) 

The  Clausius-Duhem  inequality  as  a  statement  for  second  law  of  thermodynamics  has  the  fol¬ 
lowing  local  form  for  the  representative  micro-structure 

•1/2  •k 

pVli'  -  pV  +  2  0  (6.9) 

0 

By  combining  (6.9)  and  (6.1)  and  using  (6.6)  we  have  the  inequality 
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p’(0*f|*  -  £)  +  T**j7^  -  q*'‘0i  >  0  (6.10) 

o 

which  in  teims  of  the  Helmholtz  free  energy  V|r*  defined  in  (6.8)  becomes 


-p*(V*  -t-  Ti*0*)  -I-  ^0  (6.11) 

Now  for  elastic  materials  the  constitutive  relations  for  Helmholtz  free  energy,  the  specific 
entropy  and  the  stress  tensor  can  be  expressed  in  the  following  forms 


♦  ^  *  /\*\ 

V  =V(Y«.0) 

(6.12) 

^  90* 

(6.13) 

(6.14) 

where  the  partial  derivative  with  respect  to  the  symmetric  tensor  Yy 

following  symmetric  form 

is  understood  to  have  the 

1  (iil  +  iil) 

The  constitutive  relation  for  the  heat  flux  vector  has  the  form 

q*'^  =  a;) 

(6.15) 

and  the  response  function  q**'  in  the  light  of  the  Clausius-Duhem  inequality  is  seen  to  be  res¬ 
tricted  by  the  inequality 


(6.16) 

With  the  help  of  (6.13)  and  (6.14)  the  energy  equation  (6.1)  is  reduced  to  the  following 

form 
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pV  -  q•^k  -  p’eV  =  0  (6.17) 

where  we  have  used  the  definition  (6.8)  in  order  to  calculate  E*  in  terms  of  y*  and  then  used 
the  relations  (6.13)  and  (6.14)  to  further  simplify  the  energy  equation.  It  should  be  recalled 
that  the  argument  of  different  functions  in  the  energy  equation  (6.17)  is  (0“,6^^\^,t)  and  this 
equation  is  written  for  each  and  every  representative  element  (k  =  1,2,. ..,n)  which  repeats 
itself  in  our  model  and  n  — >  <».  For  a  bi-laminate  representative  micro-structure  with  thick¬ 
ness  ^2  we  introduce  the  following  composite  quantities.  These  relations  can  be  generalized 
for  a  multi-constituent  micro-structure  without  any  difficulty  (see  definitions  5.3  and  5.4) 

,  ^ 

pg‘'^‘=f  i  P‘g'''^‘d4 

^2  0 

(6.18) 

Pg'^i  =  j  P*g*’'V^  d^ 

0 

,  ^ 

g‘''q“=f  J  g'’V‘d4 
^2  0 

(6.19) 

g''V‘=f  /g'‘V“5d5 

S2  0 

Pg'^m)^-^  J  pV''^’5”d5  (m  =  0,1,2)  (6.20) 

^2  0 

We  further  assume  that  the  variation  of  temperature  6*  across  the  micro- structure  is  a  linear 
function  of  hence 

0*(0“,03<‘‘),^,t)  =  <|)o(0“,0^^>,t)  +  ^(l)i(0“,03O‘>,t)  (6.21) 

In  order  to  derive  the  appropriate  form  of  the  energy  equation  for  the  composite  laminate,  first 
we  write  (6.17)  in  the  following  form 
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p«g*i/2r*  _  (g*i/2q*k)  =  p*g*i/2Yj*e*  (6  22) 

which  after  using  (6.21)  reduces  to 

p*g*i/V  -  (g*‘^q*‘'),k  =  p’g*^^*(<l>o  +  (6.23) 

Now  divide  (6.23)  by  ^2  and  integrate  with  respect  to  4  from  0  to  ^2.  the  result  is 

J  P*g**'Vd4  ~  {  (g**  V“).od4  “  J  (g*^V^)d^ 

^  I2  !  ^*®*'^*‘*^  ^  I2  { 

Each  term  in  the  above  equation  can  be  written  in  terms  of  the  composite  quantities  intro¬ 
duced  in  (6.18)-(6.20)  except  the  third  term  which  is  the  difference  between  the  values  of 
g*i/2  q*3  above  and  below  the  representative  element  divided  by  its  thickness  ^2»  namely 

•r  /  4*  (g*'V^)d^  =  -r  [g*''V^(0“.0^^‘‘^^^t)  -  g*''2q*3(ea  03(k)  t)j 

^2  0  ^2 

Now  we  assume  the  existence  of  the  continuous  function  h(0“,0^,t)  which  coincides  with 
q*^(0“,0^<‘'^,t)  at  0^  =  0^^\  and  further  approximate  the  right  side  of  the  above  equation  as  the 
gradient  of  this  function  multiplied  by  g^^^  in  the  0^  direction,  i.e.,  0(g^‘'^  h)/00^.  As  a  result, 
(6.24)  can  be  written  as 


pg^'^  -  (g^^q“).a  -  h)  =  pg’'^(0ono  +  <t>iT1i)  (6-25) 

In  writing  (6.25)  we  have  also  made  use  of  the  balance  of  mass  equation. 

Next  we  multiply  (6.23)  by  integrate  with  respect  to  ^  from  0  to  ^2  and  divide  it  by  %2 
to  get 
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=  y-  J  P*g*''^*^(0o  +  4<l>i)d^  (6.26) 

^2  0 

Using  integration  by  part,  the  third  term  on  the  left-hand  side  of  (6.26)  can  be  written  as 

I  4  (g*^V^)d^  =  ~  1  g*'^q*^d^  =  g^'2  (h-q3)  (6.27) 

which  in  writing  the  last  term  we  have  used  (6.19)i  and  the  definition  of  h  given  above. 
Using  this  result  together  with  the  relations  (6.18)2,  (6.19)2  and  (6.20)  we  can  write  (6.26)  in 
the  following  form 


pg^^^i  -  (g^'^q“).a  -  g'^(h-q^)  =  pg^'^(Mi  +  <l>iTq2)  (6.28) 

To  determine  the  appropriate  form  of  constraints  on  the  composite  heat  flux  vectors,  first 
we  write  the  Clausius-Duhem  inequality  (6.16)  in  the  following  form 

g*J^2  -q*ke;  ^  0 

which  by  (6.21)  reduces  to 

g*‘'2  q*“(<t)o.o  +  ^<t>i.a)  +  g*^'^  q*^<t>i  ^  0  (6.29) 

Next  we  divide  (6.29)  by  ^2  and  integrate  with  respect  to  ^  from  0  to  ^2  which  after  using 
(6.19)  can  be  written  as 


or 


g^'^  q“<t>o,a  +  g^'^  qr<l>i.a  +  g’'^  q^<t>i  ^  o 


q“<f>o.a  +  qf<!>i.a  +  q^<t>i  ^  o 


(6.30) 
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which  is  the  appropriate  form  of  Clausius-Duhem  inequality  for  the  elastic  composite  lam¬ 
inates. 

When  the  rate  of  heat  supply  or  absorption  is  zero  (r  =  r^  *  0)  the  energy  equations  for 
the  composite  reduce  to 

q“).a  +  ^  (g^^^  h)  +  p  =  0  (6.31) 

(g^'^  qr),a  +  g^^  (M^)  +  P  g^^  (<l>oili  +  =  (J  (6.32) 

For  small  deformations  of  composites  with  initially  flat  plies  the  energy  equations  (6.31)  and 
(6.32)  further  reduce  to 


r}h 

q,S  +  +  p(<l>oiio  +  =  0 

q“a  +  b  -  q^  +  p((l>oVi  +  <|)iil2)  =  0 


(6.33) 


where  obviously  no  distinction  should  be  made  between  contravariant  and  covariant  com¬ 
ponents  of  heat  flux  vectors. 


The  derivation  of  energy  equations  (6.25)-(6.28)  and  the  Clausius-Duhem  inequality 
(6.30)  for  the  composite  laminates  is  not  affected  by  the  number  of  layers  (or  constituents)  in 
the  representative  micro- structure.  The  only  necessary  modiflcation  in  the  case  of  a  multi¬ 
constituent  composite  is  the  replacement  of  ^2  by  ^  in  deflnuions  (6.18)-(6.20).  Here  ^  is 
the  thickness  of  the  representative  element  which  is  supposed  to  consist  of  n  layers.  Of 
course  ^  is  still  supposed  to  be  a  very  small  number. 


To  recapitulate,  for  a  composite  whose  micro-structure  is  composed  of  n  layers  with  a 
total  thickness  of  ^  we  have  the  following  relations  for  energy  balance  and  the  Clausius- 
Duhem  inequality 
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pr  _  g-J^[(gl'2qa)^  +  _  (gi/2  },)]  =  +  0^^^) 


pri  +  q3  -  h  -  g'^^(g‘^qi“).a  =  pCMi  +  ^1^2) 

q“<l>o.a  +  q“<l>i.a  +  q^<l>i  ^  o 
where  the  composite  quantities  are  defined  as  follows 


=  T“  J  P*g**^* 
^  0 


Pg'^^i  =  4-]  P*g**^*^ 

0 

1  ^ 

Pg^^^(in)  =  T  /  P*g*^'^*4'"d4  (m  =  0.1,2) 

0 


gl/2qi  =  J  g*l^q*‘d^ 

Sn  0 


g^'^q”  =  -r  J 
^  0 


(6.34) 


(6.35) 


7.0  CONSTITUTIVE  RELATIONS  FOR  LINEAR  THERMO-ELASTICITY 


For  a  composite  laminate  whose  micro-structure  is  composed  of  n  layers  with  different 
linear  thermo-elastic  constituents,  we  recall  the  following  constitutive  equations  for  the  stress 
tensor  entropy  and  the  heat  flux  vector  q* 


-  4)e*  (7.1) 

(P*Tl*)(a)  =  +  (P*c)(a)e*  (7.2) 

q(a)  =  -4)6  •  (7.3) 

where  C(^,  C(i),  C(5()  and  4)  (®  =  1.2,...,n)  are  constants  in  the  associated  layers.  Moreover, 
we  have  the  following  symmetries 


(7.4) 

(7.5) 


Now  we  proceed  to  calculate  the  appropriate  constitutive  relations  for  composite  stress 
vector  T‘,  composite  couple  stress  S“,  composite  entropy  Ti(n,)  (m  =  0,1,2),  and  composite  heat 
flux  vectors  q‘  and  qf.  The  contribution  of  the  first  pan  of  (7.1)  to  the  constitutive  relations 
for  T*  and  S®  (  and  consequently  and  S“j)  has  already  been  calculated  (see  section  5). 
Therefore  we  need  to  find  out  the  effect  of  the  second  pan  of  (7.1)  in  the  constitutive  rela¬ 
tions  for  T‘  and  S®.  Similar  to  what  was  dene  in  section  5  we  adopt  the  following  definitions 
for  the  weighted  averages  of  various  quantities 


j(Wij  = 


0 


J  ^‘'(P*c)(a)d^ 


(7.6) 


(7.7) 


-  82- 


4. 

=  ±  j  (7.8) 

^  0 

Now  recalling  (5.11),  (5.12),  (6.21)  and  (7.6)  we  can  write 

J  elk)  0*  gi*d^  =  "T  /  ^Ik)  g/ 

+  J  ^elk)  ^  J  (c'Pgp*  +  c^gj)  g*^'^  d^ 

0  0 

(h  ^ 

+  ~P  ]  (c‘^gf  +  C^g3)d^ 

0 

=  <l>og^^  gp  "T  J  ^pgj  -p  J  ^(1  + 

^n  Q  ^n  0  *^8 

+  <l>o  g^'^  I  gp  I  ^(1  +  -|~)c‘^d^ 

+  01  g'^  ^  gj  J  —  £■  J  ^(1  + 

Sn  0  ^8  ^  0  ^8 

=  fo  g''^  gjl-|^  I  (1  +  I  27^'''^'*^’ 

+  01  g'“  gj(  -jr  J  5('  +  +  A.J  -i  J  4*(i  +  •^)c0d5) 

^0  ^8  Sn  0  ^8 

=  ♦„  g''^  gilJ*^  +  J<"'i  +  J^O"*  +  ^  J<*')) 

+  0,  g''2  gj(J<'W  +  4-  J®"  +  +  -T  ^‘”^1 

2g  2g 


(7.9) 
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By  combining  the  results  (7.9)  and  (5.17)  we  can  obtain  the  response  function  for  the 
composite  stress  vector  T*.  As  before  if  we  disregard  the  factor  g'^  gj,  what  remains  is  the 
constitutive  relation  for  t'j  which  is  recorded  below 

x'j  =  {l(0)>jk/  +  A.  +  xi(l(l)‘pk/  +  A  i(2)ipk/)j.y,  ^ 

2g  ^  2g 

+  {I(l)ij/a  ^  A  I(2)ij/a  +  xj(l<2)ip^  +  A  l(3)iPto))K,„ 

-  ^ 

-  01 1 ^  (7.10) 

2g  2g 

Similar  steps  are  followed  to  find  the  constitutive  equations  for  S“  and  S“j.  The  contribution 
of  the  thermal  term  is 

J  i  eg,  e‘  g‘''2  g;d4  =  ♦„  g‘'^  gj  {J”)"'  +  ^ 

+  0,  g''2  g,  {JOW  +  A  jOWj  +  ^j(jO)oS  +  A  J(4)o0)j  (7.11) 

Again  combining  (7.11)  and  (5.20)  we  find  the  constitutive  relation  for  S“.  Dropping  the 
common  factor  g^^  gj  would  result  in  the  constitutive  relation  for  which  is  recorded 
below. 

S“j  =  Yy{lO)«jk/  ^  A  l(2)ojk/  ^  X-j(l(2)onfk/  +  A  l(3)a7k/)j 

2g  2g 

+  K^{l(2)aj/p  +  A  l(3)ajip  ^  X^l(3)ay/P  ^  A  i(4)ay/p)} 

-  + 4-  + 4- 

2g  ^  2g 
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-  }  (7.12) 

2g  ^  2g 

The  next  step  is  to  find  the  constitutive  relation  for  the  specific  entropies  (m  = 
0,1.2).  We  calculate  the  contribution  of  each  term  of  the  relation  (7.2)  separately.  By  (6.35)3 
we  have 


Pg*^(m)  =  7-  I  P*Tl*g*'^^"'d4 

Sn  0 

1  ^  .42 

=  -r  J  +  -^  J  (p*c)(o)e*g**^^"'d^ 

^0  ^  Sl 

Since  =  g^'^  (1  +  4^)  we  can  write  the  first  part  of  (7.13)  as 

2g 

1  ^  1/2  ^  t  A 

/  ^ikyij  8*'^^  1 1  27 

which  by  (2.31)-(2.34)  reduces  to 

4b  4b 

=  [  -^  p”(i  +  -IjWyTij 


=  +  4-  +  Aj(m+2)ipj} 

2g  2g 

The  second  part  of  (7.13),  by  (6.21)  and  (7.7),  is  written  as 

1  ^  1/2  ^  Ra 

I  (P‘c)(.)9'  g'*'^  5"<14  =  J  4”(1  +  -IlKfoHPiXp'cVid? 

(h  ^  t  A 

"  “IT  ! 

1/2  ^ 


(7.13) 


(7.14) 


(7.15) 
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=  <>0 


+  A  +  4),  {K("^1)  +  A  k("'+2)}  (7,16) 


Substituting  (7.15)  and  (7.16)  in  (7.13)  we  obtain 


pg^'^(m)  =  g^'^  +  ■—  J("^‘>y))  +  Kip(J<"«-*>‘P  +  ^  J<"«-2)iP) 


+  <|>o  g^'2  (K(m)  ^  A  K<”^^>)  +  <j>i  gl/2  (K(m+1)  ^  A  k("«-2)) 

2g  2g 


+  K<'"-^^>)  +  4>i(K<"^')  +  4-  (m  =  0,U)  (7.17) 

2g  2g 

We  start  with  (6.35)4  to  find  the  constitutive  relation  for  q*.  Substituting  from  (7.3)  in  (6.35)4 


we  get 


gl/2qi  =  y-  j  g*‘^q*'d^  =  -  J  g**'^ 


J  (1  +  ■^)(k'^e;p  +  k‘30;3)d^ 


(7.18) 


By  (6.21) 


®.p  =  <l>o.P  ^‘l^l.p 


(7.19) 


63  - 


(7.20) 
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Noting  these  results  we  simplify  (7.18).  Hence 


I  (1  + 


al/2  ^  r-K 
^  J  (1  +  ■|j)k‘3<Did^ 


=  -  g''^  J  (1  +  +  -^  J  ^(1  +  ■||)k»d| 


(7.21) 


Using  (7.8)  we  get 


gi/2  qi  =  _  gi/2  (4,^  p[L(0)iP  +  A  L«)‘P]  -h  +  A  lW] 

+  4,j[L(0)J^  AL(i)i3]} 


qi  =  -  (L(0)‘P  +  A  L<i)iP)(i,^p  -  +  A  L«)iP)<j>j  p 


_  (L(0)i3  +  A  L(l)i3)4,j 
2g 


(7.22) 


Finally  we  use  (6.35)5  derive  the  constitutive  equation  for  q“.  Similar  to  the  above 
development,  we  write 


1  ^  1  ^ 
gl/2  qO  ^  j  g*l/2q*a^  =  _  j  ^  g*l/2  k“je*d^ 
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=  -  U(i  +  +  k“30;3)d^ 

^  i  2g 

=  -  J  5(1  +  ■^)k»P((|i„,p  +  54iiji)<J5  -  -^  J  5(1  +  ■|^)k‘^<l>i<15  (7.23) 

^  0  ^8  0  *^8 

which  by  (7.8)  reduces  to 

qi“  =  -(L(‘>“P  +  A  l(2)«P)<J)  p  -  (L<2)“P  +  A  lO)“P)4>j  p  -  +  A  l(2)^)0i  (7.24) 

2g  2g  2g 

This  concludes  our  derivation  of  linear  thermo-elasdc  constitutive  relations  for  composite  lam¬ 
inates. 

For  small  deformations  of  a  composite  with  initially  flat  plies  the  foregoing  equations  are 


simplified  to  the  following  constitutive  relations: 

7ij  =  I®Uk./  +  IiSUui.«3  -  Jf  fo  -  (7.25) 

S„j  =  -  Jij'il'o  -  Ig'ifl  (7.26) 

=  +  +  (m  =  0,1,2)  (7.27) 

qi  =  -Lf ’t'o.S  -  l-4"fl.8  -  ^"’Ol  (7.28) 

qf  =  -  -  L(7»f(|,,,p  -  Ll'«(|>,  (7.29) 


The  constitutive  coefficients  ((i*!,  1(0  and  l(71  have  already  been  calculated  and  recorded  in 
equations  (5.26).  As  for  the  other  constitutive  coefficients  we  use  the  results  of  section  5. 
Comparing  definitions  (7.6)-(7.8)  with  (5.8)  and  using  the  results  (5.4)  and  (5.30)  we  can 
write 


T(k)ij  =  _1_  pk 

k+1  ,=1 


Am/ 


k+1 


(7.30) 


BASE 


(7.31) 


=  k^  li  ^ 

It  should  be  noted  that  k^J)  in  equation  (7.32)  are  the  coefficients  of  thermal  conductivity  of 
different  layers  of  the  representative  micro-structure  and  are  not  to  be  confused  with  the 
superscript  k  which  assumes  non-negative  integer  values. 

If  the  micro-structure  is  composed  of  isotropic  layers,  the  coefficients  of  thermal  stress 
c^^)  and  thermal  conductivity  k^^)  can  be  written  in  terms  of  only  one  constant  for  each  layer. 
For  such  case  we  write 


cf  =  P(oSii 

Taking  note  of  these  relations  and  relations  (7.30)  and  (7.32)  we  obtain 


(7.33) 

(7.34) 


Jf  =  6ij  L  P(,)Anv 

r=l 

^5"  =  7  I 

jf  =  y  I  PirtAm^  (7,35) 

Jf  =  7  I,  P(,)Aitv'' 

=  h  *Hr)AmT 
1=1 


=  y  ^8ij 
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(7.36) 

Consequently  the  constitutive  equations  (7.25)-(7.29)  reduce 


n  n 

Xij  =  8ij  £  -  P(r)0o)An\  +  (Uy  +  Uj.i)  L  ^(,)Anv 

1=1  1=1 

"J  (\r)Ua,a3  ~  P(r)<l>l)^^^ 

P  n 

+  Y  “j.a35ia)  iHr}^  (’^•37) 

i 

I  n  n 

Saj  =  T  4n{5aj  2:  (X(,)Uk,k  -  |3(,)(t)o)W  +  (Ua,j  +  Uj.a)  ^ 
z  r=l  ■'*'1=1 

+  J  ^^{8aj  (^r)Up.p3  - 

n 

+  (5jpUa,p3  +  Uj,a3  I  ll(r)Am,^)  (7.38) 


n 

P'n(O)  =  tP(r)Ui.i  +  (p  c)(r)0o]^*^ 

r=l 

+  Y  [P(r)Up.P3  +  (p*c)(,)<t>i]Am2  (7.39) 

1  "  • 

pn(i)  =  ^  2:  [|3(,)Ui.i  +  (p*c)(r)<l>o]Am2 

Z  r=l 

s’  Ji  +  (p*c)(,)(l>i]Am,^  (7.40) 


P^(2)  =  4  (P*C)(r)<l>o]Am3 

j  1=1 


BASE 


-90- 


+  ~  ^  [P(r)Up,p3  +  (p  c)(,)pi]Am/  (7.41) 

Qa  =  -  <l>o.a  I  k(,)AlTV  -  4  (7.42) 

r=l  z  1=1 

93  =  -  <l>i  k(,)An\  (7.43) 

9“  =  -  7  ^n<l>o.a  2  k(,)Am,2  -  1  ^  J  ^^^3  (7.44) 

Z  r=l  J  r=l 


In  relation  (7.44),  a  is  written  as  a  superscript  only  for  convenience  and  does  not  signify  the 
contravariant  position. 

The  linear  equations  of  motion  and  balance  of  energy  for  a  composite  with  initially  flat 
plies  are  derived  by  substituting  (7.25)-(7.29)  in  (2.123),  (2.124)  and  (6.33).  Using  the  results 
(5.41)  and  (5.42)  in  conjunction  with  (7.25)  and  (7.26)  we  have  the  following  equations  of 
motion  in  linear  thermo-elastic  theory 


+  l4j]pU/.ap3  -  Joj^^O.a  ~  Jaj^<t>l.a  +  +  ®j.3 

r=l 

=  iij  2:  Uj.3  L  pi^^Am,^  (7.45) 

T=1  Z  r=l 

^a]^/*^k.ta  (aj^^l.aP3  ~  Jaj^0O.a  “  ^oj^^Ua  +  —  ljji^/U|t./ 

-  Iijiu<.p3  +  Jiy^<l>o  +  +  Cj  z;  pi^^Anv 

=  7  4niij  \  +  7  ^niij.3  \  pi'^Am,^  (7.46) 

The  energy  equations  when  the  rate  of  heat  supply  or  absorption  is  zero  are  recorded  in  rela¬ 
tions  (6.33)  for  small  deformations  of  thermo-elastic  composites  with  initially  flat  plies. 
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Substituting  the  constitutive  relations  {1 2iy(J 29)  in  (6.33)  we  find  the  following  coupled 


differential  equations  for  displacement  and  temperature  fields 
(foJf  +  +  OoK*®  +  ■t'lKO'Wo 

+  (♦..Kl”  +  (|>,K®)(i>,  +  -  LS’fl.o  =  0  (7.47) 

((tj,;”  +  (JiJif  +  (liiJif’XfB  +  M>„K<‘>  +  (|>,K®)(ii„ 

+  (00*^^®  +  +  h  +  L^p*0o.p  ^3^01 

-  UKoP  -  ‘-®01.ap  -  Li3'01.o  =  0  <7-48) 

For  static  problems  in  the  absence  of  body  force  and  heat  supply,  the  foregoing  equations 
are  further  reduced  to 

U%Uk./a  +  lajipU/.paa  “  Jaj^0o.a  “  +  ^j.3  =  0  (7-49) 

lijlllUtta  +  lSi“taP3  -  41*00.0  -  4j’0l.o  +  Oj  -  'S*l"u 

-  I|>I.P3  +  J3T0O  +  41*01  =  0  (7.50) 

-  L<?0o.oP  -  l4^0i,op  -  4'*3*0i,o  =  0  (7.51) 

h  +  L|8*0o.p  +  MJ*0i.p  +  43*01  -  4KoP  -  LS0I.OP  -  1-S’0i.o  =  0  (7.52) 


Similar  to  what  was  done  previously  in  order  to  find  a  relation  between  the  director  dis¬ 
placement  and  the  gradient  of  displacement  vector,  we  enforce  the  continuity  of  the  tempera¬ 
ture  field  across  two  adjacent  micro-structures  to  derive  an  analogous  relation  between  and 
<t>i  defined  in  equation  (6.21).  In  order  that  the  temperature  field  be  continuous  on  the  com¬ 
mon  surface  between  and  (k+l)*‘  micro- structures  we  should  have 
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0*(ea  e3(k+i)  o^t)  =  6’(e“.e3<‘‘>4,t)  (7.53) 

Now  by  (6.21)  we  have 

e'(e“,e»*‘>,o,t)  =  (7,54) 

=  i|>„(6“.e^>,t)  +  ^(|),(e“,e3«,t)  (7.55) 

Substituting  from  (7.54)  and  (7.55)  in  (7.53)  we  get 

<l>„(e“,e’»*»,t)  =  <ii„(e“,e’»>,t)  +  ^(|.,(e“,e“‘>,i) 

or 

<i)i(0“,e3»>.t)  =  {(t)o(e“.e3<‘^+»),t)  -  (i),(e“,e3(w,t)}  (7.56) 

By  smoothing  assumptions  and  noting  the  smallness  of  ^  we  approximate  the  right-hand  side 
of  (7.56)  as  the  gradient  of  in  the  0^  direction.  So  we  obtain 

i|),(e»,e’,t)  =  (|i„(e“,e3,t)  (7.57) 

This  conclusion  is  used  in  various  field  equations.  In  particular,  equation  (7.49)-(7.52)  reduce 
to 

^j^/Uk./a  +  U]^U/,„p3  -  J4j^(l>0,a  ~  Jaj^‘l>0.a3  +  ^j.3  =  0  (7.58) 

^oj^t*l.oP3  “  '^aj^4^0.a  ~  •^aj^*l*0,a3  +  tJj  ~ 

“  l|jiijU/,P3  +  ”  0  (7.59) 

hj  -  L^^<l>0.aP  “  L^4>0.oP3  ~  ^^^00,03  =  0  (7.60) 

h  +  M^^^0.p3  M3^*^0.3  ~  l4V<l>0.aP  “  ^a^^0.oP3  ~  ^y4>0.a3  (7.61) 
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Eliminating  Oj  between  (7.58)  and  (7.59),  and  h  between  (7.60)  and  (7.61)  we  find  the  follow¬ 
ing  coupled  differential  equations  for  displacement  and  temperature  fields.  Since  we  are 
investigating  the  static  problems  in  the  present  derivation  the  equation  for  temperature,  i.e., 
the  equation  resulting  from  the  energy  equations  is  independent  of  the  displacement  field. 
Recalling  (5.45),  the  displacement  equation  becomes 


,la2  ~  Io^p“k.aP33 

-  jSVi  -  JiS' Vi3  +  JS  Va3  +  45Va33  =  0 


(7.62) 


and  the  temperature  equation  by  (7.60)  and  (7.61)  is 


Lf  lloaj  +  (Wi*  -  Lfi’Wo.oss  -  l^'l'0.aP33  =  0  (7.63) 

Having  determined  the  displacement  and  the  temperature  fields,  the  interlaminar  stresses  Cj 
and  heat  flux  h  can  be  determined  from  (7.59)  and  (7.61),  respectively.  The  results  are 

^Sp“k.p3  “  ^ajJ/“k.ia  ” 

*  JiiW  +  *0  -  lij’tioj  OM) 


h  -  +  L^6o.a3  “  “  Ma^a3  -  Lj3^<j>03 

=  L4V0O.ap  +  ^^<l>0,oP3  +  ~  Ma^<l>0,a3  “ 


(7.65) 
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